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This  report  is  Intended  to  provide  a  general  review  of  the 
formation,  growth,  distribution,  properties,  and  behavior  of  sea  Ice  In 
the  polar  oceans.  Special  emphasis  Is  given  to  factors  that  directly 
affect  biological  activity  beneath  the  Ice,  so  that  the  treatment  of 
some  topics  (e.g.,  mechanical  properties.  Ice  dynamics)  Is  less  complete 
than  that  of  others  (e.g.,  optical  properties.  Ice  growth).  Although 
written  primarily  for  the  nonspeclallst,  the  report  Includes  an 
extensive  bibliography  where  more  detailed  Information  can  be  obtained. 

A  slightly  modified  version  of  this  report  will  also  appear  as  Chapter 
2,  "The  Ice  environment,"  In  Sea  Ice  Biota,  edited  by  Rita  A.  Homer. 

The  material  here  Is  reprinted  with  permission  from  Sea  Ice  Biota, 
copyright  CRC  Press,  Inc.,  Boca  Raton,  Florida. 

I  am  Indebted  to  Steve  Ackley,  Don  Perovlch,  and  Ed  Vaddlngton  for 
technically  reviewing  the  material  In  this  report.  Special  thanks  go  to 
Rita  Homer  for  her  patience  and  many  helpful  suggestions  during 
preparation  of  the  manuscript.  Finally,  1  would  like  to  acknowledge  the 
contribution  made  by  the  Office  of  Naval  Research  whose  long-term 
commitment  to  sea  Ice  research  produced  much  of  the  Information 
summarized  here,  and  whose  support  under  Contract  N00014-84-C-0111  made 
possible  this  review. 
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1.  EXTENT  AND  MORPHOLOGY  OF  THE  POLAR  ICE  PACK 


Sea  ice  Interacts  thermodynaalcally  with  both  the  atmosphere  and 
ocean,  giving  rise  to  a  climate  over  the  polar  oceans  that  is  more  char¬ 
acteristic  of  the  continental  ice  sheets  than  of  a  marine  environment. 

By  area,  sea  ice  accounts  for  roughly  two-thirds  of  the  earth's  perma¬ 
nent  ice  cover,  but  only  about  0.1%  of  its  volume.  Thus,  unlike  the 
polar  ice  caps,  sea  ice  is  only  a  thin  veneer  whose  thickness  and  extent 
can  undergo  large  changes  in  response  to  small  changes  in  climate  or 
oceanic  heat  transport.  Not  only  does  a  sea  ice  cover  greatly  reduce 
heat  exchange  between  the  surface  and  the  atmospheric  boundary  layer, 
but  it  also  affects  the  temperature  and  salinity  structure  of  the  under¬ 
lying  ocean.  In  particular,  growth  and  melting  of  the  ice  causes 
continual  salinity  changes  at  the  top  of  the  mixed  layer,  while  varia¬ 
tions  in  ice  thickness  and  concentration  complicate  the  input  of  solar 
energy  to  the  water.  Because  of  its  impact  on  the  ocean  and  atmosphere 
and  because  of  a  positive  feedback  between  ice  extent  and  surface 
albedo,  sea  ice  is  generally  considered  to  be  an  important  element  in 
the  earth's  climate  system. 

A.  Geographical  Distribution 

At  maximum  extent  sea  ice  covers  some  8%  of  the  Southern  Hemisphere 
and  5%  of  the  Northern  Hemisphere.  The  nature  of  the  ice  pack  in  the 
two  hemispheres  is,  however,  quite  different  primarily  because  of 
differences  in  the  distribution  of  the  polar  land  masses.  In  the  south 
the  Antarctic  Continent  occupies  most  of  the  area  poleward  of  70°S,  so 
that  sea  ice  forms  a  seasonally  varying  annulus  in  the  region  between 
55°  and  70°S  (Fig.  1).  Meridional  variations  are  relatively  small,  the 
greatest  northward  extent  occurring  in  the  Weddell  Sea  sector.  Most  of 

C  O 

the  antarctic  ice  pack  is  seasonal,  covering  20 '10  km  at  maximum 
extent  and  only  4 ‘10  km  at  the  end  of  summer  (Zwally  et  al.,  1983). 

In  the  Northern  Hemisphere  the  area  poleward  of  70°N  is  in  large 
part  a  landlocked  ocean  basin  which  is  ice  covered  during  most  of  the 
year.  Seasonal  ice  forms  in  many  of  the  peripheral  seas  surrounding  the 


Arctic  Basin  (Fig.  2).  At  minimum  extent  (August)  the  Ice  covers  about 
6  2 

7 ‘10  km  ,  almost  double  that  in  the  Antarctic,  while  at  maximum  extent 
(March)  the  area  expands  to  14*10^  km^  (Walsh  and  Johnson,  1979). 


kilometers 
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The  extent  of  sea  Ice  In  the  Southern  Hemisphere  with 
greater  than  1/8  concentration  (after  CIA,  1978). 
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greater  than  1/8  concentration  (after  CIA,  1978) 


Variations  in  oceanic  circulation  and  the  complex  configuration  of  the 


land  masses  cause  meridional  variations  in  winter  ice  extent  of  up  to  35 
degrees  of  latitude.  Cold  air  flowing  over  the  Sea  of  Okhotsk  from  the 
Asian  Continent  results  in  the  formation  of  a  seasonal  ice  cover  which 
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typically  extends  to  the  north  coast  of  Japan  (44°-45°N) .  At  the  other 
extreme  Is  the  area  west  of  Spitsbergen  (0°,  80^N)  where  warm  Atlantic 
water  flowing  Into  the  Arctic  Basin  and  the  southward  drift  of  the  Ice 
combine  to  maintain  the  Ice  edge  at  a  nearly  constant  position  through¬ 
out  the  year.  The  economic,  military,  and  scientific  Importance  of  the 
region,  together  with  Its  relative  accessibility,  have  stimulated  exten¬ 
sive  research  efforts  and  we  now  have  a  good  general  picture  of  the  Ice- 
ocean-  atmosphere  system  In  the  Arctic.  Such  Is  not  the  case  In  the 
Antarctic.  Although  some  satellite-derived  data  are  available  regarding 
the  drift,  extent,  and  physical  characteristics  of  antarctic  pack  Ice 
(e.g.,  Ackley,  1979),  conditions  In  most  Ice  covered  portions  of  the 
Southern  Ocean  are  poorly  known. 


Sea  Ice,  driven  principally  by  winds  and  currents.  Is  In  almost 
continual  motion.  As  will  be  discussed  later,  these  motions  affect  not 
only  the  extent  of  the  ice,  but  also  its  physical  properties  and  the  way 
In  which  It  Interacts  with  the  atmosphere  and  ocean.  Ice  motion  In  the 
Arctic  Basin  has  been  well  documented  from  the  trajectories  of  manned 
drifting  stations  and  automatic  data  buoys.  Most  references  describe 
the  two  major  drift  features  Identified  by  Gordienko  (1958):  (1)  the 

Transpolar  Drift  Stream  which  transports  Ice  from  the  western  side  of 
the  basin  across  the  pole  and  through  Fram  Strait  between  Greenland  and 
Spitsbergen,  and  (11)  a  clockwise  gyre  in  the  Beaufort  Sea  region.  How¬ 
ever,  a  recent  analysis  of  the  mean  drift  field  by  Colony  and  Thorndike 
(1983)  shows  a  somewhat  different  pattern  with  antlcyclonlc  motions 
throughout  the  Ameraslan  Basin  and  southward  drift  toward  Fram  Strait  in 
the  Eurasian  Basin.  Average  displacements  are  approximately  7  km/day. 
Ice  movement  In  the  Greenland  and  Bering  seas  tends  to  be  southerly  with 
average  velocities  on  the  order  of  15  km/day  (Muench  and  Ahlnas,  1976). 
While  some  Ice  Is  lost  from  the  Arctic  Basin  through  the  Bering  Strait, 
this  Is  small  compared  to  the  export  through  Fram  Strait.  On  an  annual 
basis  roughly  10%  of  the  Ice  In  the  basin  drifts  Into  the  Greenland  Sea 
where  it  eventually  melts,  providing  a  significant  heat  gain  for  the 


basin.  Ice  movements  In  the  Antarctic  are  much  simpler,  being  generally 
clockwise  around  the  continent.  A  major  drift  feature  exists  In  the 
Weddell  Sea  In  the  form  of  a  cyclonic  gyre  which  advects  Ice  along  the 
Antarctic  Peninsula  and  allows  multiyear  Ice  to  develop  In  the  region. 

Diverging  motions  In  the  Ice  cause  it  to  break  apart  and  expose  the 
\inderlylng  water.  These  cracks,  called  leads,  tend  to  be  narrow  (tens 
to  hvindreds  of  meters  In  width) ,  but  may  extend  for  many  kilometers .  In 
some  situations  leads  may  reach  several  kilometers  in  width  and  are 
termed  polynyas.  The  steep  temperature  gradients  which  exist  over 
winter  leads  result  In  large  turbulent  heat  losses  to  the  atmosphere  and 
rapid  Ice  production.  During  the  summer,  leads  act  as  moving  windows 
which  admit  substantial  amounts  of  solar  energy  to  the  upper  ocean. 

Leads  are  common  In  most  regions.  In  the  Arctic  Basin  winter  leads 
occur  with  an  average  spacing  of  about  300  m.  Leads  tend  to  freeze  over 
quickly  and  open  water  rarely  makes  up  more  than  0.5%  of  the  winter  Ice 
cover.  Young  ice  (h  <  1  m)  in  refrozen  leads  typically  accounts  for  5- 
10%  of  the  area.  During  the  summer  decay  cycle  the  area  of  open  water 
increases  greatly,  reaching  maxlmvua  values  of  10-20%  in  the  Central 
Arctic. 

In  marginal  seas  where  the  Ice  edge  Is  not  confined  by  land,  ice 
concentration  patterns  are  complex  and  can  vary  rapidly  with  time.  Over 
the  Interior  of  the  Bering  Sea  Ice  pack,  for  example.  Cavalier!  et  al. 
(1983)  found  winter  concentrations  ranging  between  85%  and  100%  with 
cold  northeast  winds,  but  these  values  dropped  to  as  low  as  45% 
following  the  onset  of  warm  southerly  winds.  Ice  concentrations  In  the 
Antarctic  are  also  fairly  low,  averaging  60-80%  throughout  most  of  the 
year  {Zwally  et  al.,  1983).  The  large  amount  of  open  water  suggests 
that  antarctic  Ice  undergoes  substantial  divergence.  In  fact,  analysis 
of  satellite  Imagery  Indicates  that  the  antarctic  Ice  pack  Is  everywhere 
divergent  In  Its  mean  motion  and  that  northward  advectlon  of  the  Ice 
must  result  In  the  continual  production  of  new  leads  and  polynyas  (Budd, 
1975;  Ackley  and  Keliher,  1976). 

Convergence  and  shearing  motions  cause  the  Ice  to  fracture  and  pile 
up  Into  a  variety  of  deformation  features,  the  most  notable  being  pres- 


sure  ridges.  Pressure  ridges  form  linear  features  up  to  several  kilo¬ 
meters  In  length  which  consist  of  a  keel  beneath  the  Ice  and  a  sail  on 
the  surface  of  the  Ice.  Keels  are  typically  about  five  times  thicker 
than  sails.  Most  ridges  In  the  Central  Arctic  appear  to  be  built  out  of 
Ice  that  was  originally  less  than  1  m  thick  (Herbert,  1970).  This  may 
not  be  true  In  all  areas,  however,  as  detailed  morphological  studies 
north  of  Prudhoe  Bay,  Alaska,  Indicated  that  nearly  50%  of  the  ridges 
there  were  composed  of  blocks  between  1.0  and  1.6m  thick  (Tucker 
and  Govoni,  1981).  Ridge  frequency  In  the  Arctic  Basin  varies  with 
season  and  location,  but  generally  ranges  from  two  to  six  per  kilometer. 
Keel  depths  are  usually  less  than  20  m.  It  has  been  estimated  (Wadhams, 
1980)  that  as  much  as  40%  of  the  total  mass  In  the  arctic  Ice  pack  may 
be  contained  In  deformation  structures.  In  contrast  to  the  arctic  situ¬ 
ation,  ridges  are  generally  believed  to  be  uncommon  In  most  parts  of  the 
Antarctic  because  of  the  divergent  nature  of  the  velocity  field  there. 

A  notable  exception  Is  the  Weddell  Sea  where  ridged  Ice  may  accotint  for 
more  than  15%  of  the  area  (Ackley,  1979) . 

Without  continued  deformations  the  polar  Ice  pack  would  assume  a 
much  more  uniform  character.  In  areas  of  seasonal  Ice  the  thickness 
would  Increase  smoothly  Inward  from  the  edge,  while  In  areas  of  peren¬ 
nial  Ice  the  thickness  within  a  particular  climatic  regime  would  be 
essentially  constant.  The  reality  Is  that  a  typical  sea  Ice  cover 
contains  a  great  range  of  thicknesses,  all  temporarily  coexisting  under 
similar  thermal  forcing.  Thicknesses  tend  to  be  distributed  chaotically 
with  large  variations  occurring  on  horizontal  scales  of  only  a  few 
meters.  Ice  motion  enhances  the  transfer  of  momentum  from  the  atmo¬ 
sphere  to  the  ocean  and  causes  thermodynamic  Interactions  between  the 
air  and  water  to  be  much  more  vigorous  than  would  otherwise  be  the  case. 
Mechanical  breakup  and  advectlon  also  act  to  accelerate  the  decay  and 
retreat  of  the  Ice  cover,  allowing  larger  portions  of  the  polar  oceans 
to  become  Ice  free  each  summer. 

C.  Ice  Zones 

Because  of  differences  In  the  properties  of  the  Ice,  In  the 
dominant  Ice-ocean  Interaction  processes,  and  In  the  nature  of  the 


scientific  problems  encountered,  the  polar  Ice  pack  Is  frequently 
subdivided  Into  different  zones.  The  most  basic  distinction  Is  between 
regions  that  are  perennially  Ice  covered  and  those  with  seasonal  Ice. 
Other  ice  regimes  frequently  mentioned  In  the  literature  Include  the 
fast  Lee  zone,  the  shear  zone,  and  the  marginal  ice  zone. 

1 .  Perennial  Ice  Zone 

Areas  of  year-round  Ice  are  found  In  both  the  Arctic  and  Antarctic; 
however,  the  characteristics  of  this  Ice  differ  substantially  In  the  two 
hemispheres.  The  perennial  ice  zone  In  the  Arctic,  located  largely 
poleward  of  75-80°N  (Fig.  2),  contains  about  two-thirds  of  all  the 
multiyear  sea  Ice  found  In  the  World  Ocean.  Thickness  averages  3-4  m  In 
the  central  part  of  the  basin,  but  Increases  significantly  north  of 
Greenland  and  Ellesmere  Island  as  the  motion  of  the  Ice  toward  the  land 
causes  extensive  ridging.  The  age  of  this  Ice  Is  uncertain,  but  may 
reach  tens  of  years  In  the  Beaufort  sector  and  at  least  several  years  on 
the  Soviet  side  of  the  basin. 

Multiyear  Ice  In  the  Antarctic  Is  restricted  to  the  Weddell  Sea  and 
parts  of  the  Ross,  Amundsen,  and  Bellingshausen  seas  (Fig.  1).  The  age 
of  this  Ice  Is  probably  much  less  than  In  the  Arctic.  What  is  known  of 
drift  patterns  In  the  Weddell  Sea  suggests  that  the  residence  time  of 
the  Ice  Is  not  more  than  2-3  years.  Nevertheless,  Ice  thicknesses  of  up 
to  5  m  have  been  observed  In  undeformed  Ice  there  (Weeks  and  Ackley, 
1982).  The  structure  of  this  Ice,  however,  was  quite  different  from 
that  In  the  Central  Arctic.  Some  60-70%  of  the  thickness  was  made  up  of 
layers  of  small  Ice  crystals  a  few  millimeters  across  rather  than  the 
organized  columnar  structure  commonly  observed  In  most  arctic  Ice. 

These  so-called  "frazil"  crystals  could  originate  In  the  water  column 
beneath  the  Ice  or  In  nearby  leads,  but  the  exact  mechanisms  Involved 
have  not  yet  been  positively  Identified.  In  addition  to  the  more 
complex  Internal  structure,  weak  surface  melting  during  the  antarctic 
summer  allows  the  salinity  of  the  ice  near  the  surface  to  remain  high. 
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As  a  result,  microwave  techniques  that  rely  on  differences  In  near¬ 
surface  salinity  to  distinguish  between  first-year  and  multiyear  Ice  are 
much  less  successful  In  the  Antarctic  than  In  the  Arctic  (see  Section 
4E). 

2.  Seasonal  Sea  Ice  Zone 

6  2 

The  seasonal  sea  ice  zone  (SSIZ)  covers  about  25*10  km  ,  roughly 
7%  of  the  World  Ocean.  The  area  covered  by  the  SSIZ  In  the  Arctic  Is 
about  the  same  as  that  of  the  perennial  Ice.  While  undeformed  Ice  In 
the  arctic  SSIZ  Is  generally  less  than  2  m  thick,  large  amounts  of 
ridged  Ice  are  present  so  that  average  thicknesses  can  be  much  greater, 
particularly  near  the  coasts.  As  pointed  out  previously,  80%  of  the 
highly  mobile  antarctic  Ice  pack  Is  seasonal.  It  has  been  reported 
(Paige,  1966)  that  seasonal  Ice  In  the  Antarctic  can  grow  to  thicknesses 
at  least  50%  greater  than  In  the  Arctic.  This  Ice  has  frequently  been 
observed  to  possess  not  only  layers  of  frazil  crystals,  but  also  a  layer 
of  Infiltration  Ice  formed  by  flooding  of  the  Ice  surface  In  response  to 
snow  loading  In  the  early  stages  of  growth.  The  generality  of  these 
observations,  however,  has  not  been  established,  and  It  Is  not  known 
whether  they  apply  throughout  the  southern  SSIZ  or  only  to  certain 
regions  near  the  continent.  The  relative  Importance  of  various  Ice 
decay  processes  also  differs  In  the  two  hemispheres.  Surface  melting, 
for  example.  Is  a  major  factor  In  the  summer  disappearance  of  arctic 
Ice.  Antarctic  Ice,  on  the  other  hand.  Is  characterized  by  minimal 
surface  melting  and  an  absence  of  melt  ponds,  probably  owing  to  lower 
relative  humidity  In  the  atmosphere  (Andreas  and  Ackley,  1982).  This 
means  that  heat  from  the  water  must  be  largely  responsible  for  seasonal 
Ice  retreat  In  the  Southern  Ocean.  In  addition,  higher  biological 
activity  In  antarctic  Ice  enhances  Internal  melting,  weakening  the  Ice 
and  accelerating  the  breakup. 

Ice -ocean- atmosphere  Interactions  In  the  shelf  regions  of  the  SSIZ 
appear  to  have  an  Important  effect  on  the  large-scale  structure  and 
circulation  of  the  World  Ocean  (e.g..  Gill,  1973).  It  has  long  been 
suspected  that  salt  fractionation  during  freezing  leads  to  the  formation 


of  cold,  dense  water  on  the  shelves.  Such  water  appears  to  be  the  pre¬ 
cursor  to  the  Antarctic  Bottom  Vater  which  fills  much  of  the  deep  World 
Ocean.  Few  details  are  known  regarding  production  sites,  rates,  or 
mechanisms  for  the  formation  of  this  water.  Offshore  winds  produce  a 
number  of  large  recurrent  polynyas  arotmd  the  Antarctic  Continent,  espe¬ 
cially  In  the  Ross,  Weddell,  and  Amundsen  seas  {Streten,  1973;  Carsey, 
1980) .  Rapid  Ice  growth  In  such  polynyas  and  the  resulting  Increase  In 
density  leads  to  vertical  hallne  convection,  possibly  all  the  way  to  the 
sea  floor  In  some  areas.  It  Is  probable  that  much  of  the  Antarctic 
Bottom  Water  originates  In  this  way.  Similar  processes  on  certain  of 
the  arctic  shelves  also  appear  to  produce  dense  water  which  Is  critical 
In  maintaining  the  themohallne  structure  of  the  Arctic  Ocean  (HcPhee, 
1980). 

3 .  Fast  Ice  Zone 

Shorefast  Ice  occurs  along  most  coasts  In  the  Arctic  and  Antarctic. 
It  forms  early  In  the  winter  In  shallow  water  where  the  water  column  can 
cool  rapidly  to  the  freezing  point.  It  Is  anchored  to  the  bottom  at  the 
shoreline  and  frequently  at  depths  of  up  to  20  m  by  grounded  pressure 
ridge  keels  {Kovacs  and  Hellor,  1974).  Although  rigidly  anchored,  fast 
Ice  may  move  tens  of  meters  In  response  to  theinnal  and  mechanical 
stresses  {Tucker  et  al. ,  1980) .  Such  motions  are  Important  because  of 
their  potential  Impact  on  offshore  structures.  Salinities  beneath  Ice 
In  lagoons  and  basins  where  circulation  Is  restricted  can  reach  values 
In  excess  of  100°/oo.  The  extent  of  fast  Ice  Is  determined  primarily  by 
sea  bottom  and  shoreline  topography  and  Is  highly  variable.  Near  Point 
Barrow,  for  example,  the  zone  Is  about  15  km  In  width  while  at  nearby 
Harrison  Bay  It  extends  some  60  km  offshore.  Fast  Ice  development  Is 
most  extensive  In  areas  where  the  shelf  slope  Is  gentle  (e.g.,  Mackenzie 
Delta  and  Laptev  Sea) .  Also  usually  classified  as  landfast  Is  much  of 
the  Ice  that  fills  constricted  bays  and  channels  In  regions  like  the 
Canadian  Archipelago.  A  review  of  landfast  Ice  and  biological  charac¬ 
teristics  of  the  underlying  water  has  recently  been  given  by  Newbury 
(1983). 


4 .  Shear  Zone 

Aircraft  and  submarine  sonar  observations  Indicate  the  existence  of 
a  zone  of  highly  deformed  Ice  In  coastal  regions  along  the  Alaskan  North 
Slope  and  Canadian  Archipelago  (Tucker  et  al.,  1979;  Weeks  and  Ackley, 
1982).  The  origin  of  this  so-called  shear  zone  Is  landward  movement  of 
the  Ice  pack,  causing  the  Ice  to  converge,  shear,  and  ridge.  The  fre¬ 
quency  of  ridges  In  the  shear  zone  Is  several  times  greater  than  In  the 
central  pack,  but  keel  depths  are  not  appreciably  larger.  The  size  of 
the  shear  zone  depends  on  definition.  Wadhams  (1980)  defines  the  shear 
zone  as  the  region  where  Ice  characteristics  differ  significantly  from 
those  In  the  Central  Arctic.  His  analysis  of  submarine  data  suggests 
that  the  shear  zone  extends  some  400  km  off  the  coast  of  Ellesmere 
Island  where  Ice  motion  Is  normal  to  the  coast,  and  about  160  km  off 
Alaska  where  the  motion  Is  more  parallel  to  the  land.  Microwave  Imagery 
(Campbell  et  al. ,  1976)  shows  that  the  landward  part  of  this  zone  Is 
made  up  largely  of  first-year  Ice  and  that,  as  might  be  expected,  the 
outer  part  contains  a  substantial  amount  of  multiyear  Ice.  The  shear 
zone  Is  clearly  a  boundary  layer  between  the  fast  Ice  and  central  pack 
and,  as  such,  understanding  Its  behavior  Is  of  particular  Importance  In 
efforts  to  develop  large-scale  dynamic  models  of  sea  Ice. 

5.  Marginal  Ice  Zone 

One  of  the  most  complex  and  scientifically  Interesting  regions,  the 
marginal  ice  zone  (MIZ) ,  Is  located  near  the  boundary  between  the  Ice 
and  open  ocean.  Here,  penetration  of  surface  waves  Into  the  pack  breaks 
the  Ice  Into  numerous  small  floes  whose  average  size  Increases  rapidly 
with  distance  from  the  edge.  In  regions  like  the  Bering  Sea,  where 
floes  are  usually  thin  (<  0.5  m) ,  Ice  In  the  outer  5-15  km  Is  often 
heavily  rafted  and  many  reach  thicknesses  of  up  to  5  m  (Bauer  and 
Martin,  1980).  The  MIZ  Is  usually  taken  to  Include  that  part  of  the  Ice 
that  Is  In  some  way  affected  by  the  presence  of  the  open  ocean,  a  region 
some  150-200  km  In  width.  The  outer  part  of  the  MIZ  Is  characterized  by 
large  horizontal  gradients  In  the  properties  of  the  Ice,  ocean,  and 


atmosphere.  Advectlon  of  air  across  the  ice  edge  can  produce  large 
changes  In  turbulent  heat  transfer,  surface  stress,  cloudiness,  and  the 
radiative  fluxes.  Comparable  changes  also  occur  In  the  mechanical 
properties  of  the  Ice,  surface  rou^ness,  and  solar  Input  to  the  ocean. 
Horizontal  variations  In  density  structure  and  surface  stress  can  give 
rise  to  a  variety  of  mesoscale  phenomena  In  the  ocean  (eddies,  jets, 
upwelllng)  which  Impact  biological  activity  and  acoustical  properties 
across  the  MIZ.  During  the  svimmer,  floe  breakup  and  Increasing 
stratification  beneath  the  Ice  affect  the  response  of  the  Ice  to  winds 
and  currents  and  the  rate  at  which  heat  Is  transferred  from  the  water  to 
the  Ice.  Biological  activity  In  the  MIZ  Is  high  (Alexander,  1980). 
Substantial  research  efforts  are  currently  under  way  In  the  Bering  and 
Greenland  seas  to  Investigate  not  only  the  properties  of  the  Ice  cover, 
but  also  oceanographic,  atmospheric,  biological,  and  acoustical 
characteristics  of  the  MIZ  (NcPhee,  1983;  MIZEX  West  Study  Group,  1983). 
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2.  FORMATION  AND  GROWTH  OF  SEA  ICE 


The  presence  of  salt  causes  seawater  to  respond  somewhat  differ¬ 
ently  to  atmospheric  cooling  than  freshwater.  A  density  vs.  temperature 
plot  of  freshwater  (Fig.  3a)  shows  the  existence  of  a  density  maximum  at 
3.98°C.  Surface  cooling  in  a  freshwater  body  whose  temperature  is  above 
this  value  causes  a  density  increase,  generating  vertical  convection 
that  continues  until  all  the  water  reaches  the  temperature  of  maximum 
density  (T^) .  Further  cooling  then  produces  a  density  decrease  in  the 
surface  layer,  allowing  the  water  to  become  thermally  stratified.  Since 
conduction  rather  than  convection  is  then  the  dominant  heat  transport 
mechanism  in  the  water,  there  is  a  sharp  decrease  in  the  upward  flux  of 
heat  to  the  surface,  and  cooling  proceeds  rapidly.  Ice  can  thus  form  in 
freshwater  when  most  of  the  water  column  is  significantly  above  the 
freezing  point. 
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Fig.  3a.  Density  vs.  temperature  diagram  for  freshwater  and  ice 
at  atmospheric  pressure. 


Addition  of  salt  to  the  water  depresses  the  freezing  point  of  the 
solution  (T^)  according  to  the  relation  (Neumann  and  Pierson,  1966) 

Tf  -  -0.003  -  0.0527S„  -  0.000048^^,  (1) 

where  temperature  Is  in  °C  and  the  salinity  of  the  water  (S^)  is  in 
°/oo.  For  many  purposes  Eq.  (1)  can  be  approximated  simply  by  T£  - 
•0.0558^.  8alt  also  reduces  but  at  a  more  rapid  rate  than  T£ 

(Fig.  3b).  If  8^  >  24.7  °/oo,  Tj^  -  T£.  As  this  is  the  case  in  normal 
seawater,  surface  cooling  results  in  a  density  increase  and  vertical 
mixing  which  continues  until  the  water  reaches  T£.  Unlike  freshwater 
lakes,  however,  density  gradients  in  the  upper  ocean  usually  limit  the 
depth  (Zj.)  to  which  the  water  must  be  cooled  before  freezing  can 
commence.  8oviet  data  (Doronin  and  Kheisin,  1975)  Indicate  that  z^  is 
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Fig.  3b.  Effect  of  salinity  on  the  temperature  of  maximum  density 
(dashed  curve)  and  the  freezing  point  temperature  (solid 
curve) . 


usually  In  the  10-40  m  range,  although  values  in  excess  of  70  m  have 
been  observed.  The  density  structure  of  the  ocean  Is  thus  a  major 
factor  In  determining  the  onset  of  freezing.  Observations  (.Doronin  and 
Kheisin,  1975)  show  a  delay  of  as  much  as  two  months  in  the  date  of 
initial  ice  formation  in  regions  where  z^  •  50  m  compared  to  those  where 
Zg  -  10  m. 

B.  Initial  Ice  Formation 

Once  the  upper  ocean  reaches  the  freezing  point,  additional  heat 
loss  produces  slight  supercooling  of  the  water  and  ice  formation  begins. 
The  amount  of  supercooling  necessary  to  initiate  ice  growth  is  usually 
small,  probably  on  the  order  of  a  few  hundredths  of  a  degree  Celsius, 
although  observations  in  some  antarctic  polynyas  and  near  Greenland  have 
shown  a  supercooling  of  0.2-0.4°C  down  to  depths  of  tens  of  meters 
{Doronin  and  Kheisin,  1975).  Initial  ice  formation  occurs  at  or  near 
the  surface  of  the  water  in  the  form  of  small  platelets  and  needles 
which  are  termed  frazil  ice.  Frazil  crystals  do  not  generally  exceed 
3.4  mm  in  diameter.  Continued  freezing  results  in  the  production  of 
grease  ice,  a  soupy  mixture  of  unconsolidated  frazil  crystals.  Because 
grease  ice  suppresses  capillary  wave  action,  it  appears  as  a  smooth, 
grayish  layer  whose  albedo  is  only  slightly  above  that  of  open  water. 
Under  quiescent  conditions  the  frazil  crystals  quickly  freeze  together 
to  form  a  solid,  continuous  ice  cover  with  a  thickness  between  1  and 
10  cm.  More  typically,  wind  waves  induce  turbulence  in  the  water  and 
abrasion  between  the  crystals  which  Inhibit  development  of  a  solid 
cover.  Wind  and  wave  action  advect  frazil  crystals  downwind  where 
accumulations  up  to  1  m  thick  can  form  in  front  of  obstacles  (e.g. ,  floe 
edges).  Laboratory  experiments  (Martin  and  Kauffman,  1981)  indicate 
that  when  the  ice  volume  is  low,  grease  ice  behaves  as  a  fluid  whose 
viscosity  increases  as  the  proportion  of  ice  Increases.  Once  the  ice 
fraction  exceeds  30-40%,  sintering  and  regelation  create  sufficient 
bonding  between  individual  crystals  that  their  mobility  is  drastically 
reduced  and  the  transition  to  a  solid  cover  begins.  In  the  presence  of 


a  wave  field,  this  transition  is  usually  marked  by  the  formation  of 
pancakes,  rounded  masses  of  semiconsolidated  slush  0.3 *3.0  m  in 
diameter.  Pancakes  often  display  Irregular  ice  rims  around  their 
perimeters,  built  up  as  a  result  of  continual  bumping  into  one  another. 
With  time,  the  pancakes  consolidate  and  are  welded  together  into  a 
composite  ice  sheet  by  the  freezing  of  grease  ice  that  forms  between 
them.  A  detailed  review  of  Initial  freezing  processes,  together  with 
photographs  of  different  stages  of  ice  development,  appears  in  Weeks  and 
Ackley  (1982). 

C .  Young  Ice  Growth 

In  the  early  stages  of  growth  (grease  ice,  pancake  ice)  the  surface 
of  the  ice  is  wet  and  its  temperature  stays  close  to  the  freezing  point. 
With  air  temperatures  often  in  the  -10  to  -20°C  range,  turbulent  and 
longwave  radiation  losses  are  large  and  ice  formation  rapid.  However, 
once  the  ice  becomes  consolidated,  free  water  is  no  longer  available  at 
the  surface  and  its  temperature  begin^^  to  approach  that  of  the 
atmosphere.  Additional  growth  then  occurs  by  accretion  on  the  underside 
of  the  sheet,  the  rate  depending  on  how  rapidly  heat  can  be  conducted 
from  the  ice-water  Interface  toward  the  surface,  i.e.,  on  the  thermal 
conductivity  and  temperature  gradient  in  the  ice. 

While  winds,  cloudiness,  relative  humidity,  snowfall,  and  oceanic 
heat  flux  all  affect  the  rate  of  ice  growth,  field  measurements  suggest 
that  reasonable  growth  estimates  can  be  obtained  from  air  temperature 
data  alone.  Figure  4  shows  ice  thicknesses  observed  by  Anderson  (1961) 
near  Thule,  Greenland,  plotted  against  the  cumulative  number  of 
freezing- degree  days  6,  where 

rt 

e  =  (Tf  -  T^)dt  ,  (2) 

t  is  time,  and  is  the  air  temperature.  Anderson  found  that  in  the 
thickness  range  10  <  h  <  80  cm  his  data  could  be  fit  by  a  simple  power 
law, 

h^  +  5.1  h  -  6.70  .  (3) 
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CUMULATIVE  FREEZING  DEGREE  (^C)  DAYS 

Fig.  4.  Relationship  between  the  thickness  of  young  sea  ice  with 
minimal  snow  cover  and  the  cumulative  number  of  freezing- 
degree  days  (after  Anderson,  1961). 

As  the  ice  becomes  thicker,  the  close  relationship  between  6  and  h 
begins  to  break  down.  Figure  4,  for  example,  shows  that  when  6  -  3000 
degree -days,  the  observed  h  varies  between  90  and  140  cm.  Factors 
contributing  to  this  breakdown  are  the  Increasing  thermal  mass  of  the 
ice,  the  decreasing  magnitude  of  the  heat  flux  in  the  ice  relative  to 
that  in  the  ocean,  and  increasing  snow  cover.  The  term  young  ice  has 
been  somewhat  narrowly  defined  in  the  Glossary  of  Snow  and  Ice 
(Armstrong  et  al.,  1973)  as  ice  in  the  10-30  cm  category.  Here  the  term 
will  be  used  in  a  broader  sense  to  refer  to  ice  whose  temperature 
profile  is  nearly  linear  and  whose  growth  can  usually  be  described  by 
some  simple  function  of  6,  i.e.,  h  <  80-100  cm. 
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Alternate  expressions  for  h(8)  have  been  derived  by  investigators 
working  in  other  regions  (see  Bilello,  1961).  Implicit  in  such  formu¬ 
lations  is  the  idea  that  the  exact  form  of  T^(t)  is  unimportant  in 
determining  h;  rather  it  is  only  the  area  under  the  curve  that  needs  to 
be  known.  This  is  equivalent  to  saying  that  the  ice  would,  for  example 
attain  the  same  thickness  in  10  days  at  -20^C  as  it  would  in  20  days  at 
-10°C  (assuming  that  Tj  -  0°C).  Differences  in  growth  predicted  by  the 
various  empirical  formulas  presumably  reflect  differences  in  environ¬ 
mental  parameters  (snowfall,  heat  content  of  the  underlying  water 
column,  etc.)  which  are  not  explicitly  treated  in  the  formulas.  As  a 
result  such  equations  must  be  considered  at  least  partially  site- 
specific  and  their  utilization  approached  with  caution.  Applying  an 
equation  obtained  from  data  in  the  Central  Arctic  to  coastal  conditions 
or  to  conditions  in  the  Southern  Ocean  could  lead  to  substantial  errors 
To  establish  a  more  general  relationship,  it  is  necessary  to 
consider  the  physics  governing  the  situation.  The  usual  approach  is  to 
formulate  a  model  using  energy  balance  equations  at  the  top  and  bottom 
of  the  ice  (see  Section  5B)  which  are  coupled  by  an  equation  describing 
the  rate  of  heat  transport  between  the  two  boundaries.  The  problem  is 
greatly  simplified  in  the  case  of  young  ice  because  the  conductive  flux 
is  directly  proportional  to  the  temperature  difference  between  the 
boundaries.  It  can  easily  be  shown  (e.g.,  Maykut,  in  press)  that  the 
form  of  Anderson's  equation  follows  from  the  assumption  of  a  linear 
temperature  gradient  in  the  ice  and  negligible  heat  flux  in  the  water. 
The  theoretical  expression  can  be  written 

2k  2k.  2ki 
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where  h^  is  the  initial  thickness  of  the  ice,  k^  is  the  average  thermal 
conductivity  of  the  ice,  k^  Is  the  thermal  conductivity  of  the  snow, 
is  the  density  of  the  ice,  L  is  the  latent  heat  of  fusion,  is  a  bulk 
transfer  coefficient  parameterizing  turbulent  heat  transfer  between  the 


ice  and  atmosphere,  and  is  the  snow  depth  which  is  assumed  to  remain 
constant  over  the  period  of  integration.  The  net  radiation  balance  is 
assumed  to  be  zero.  When  h^  -  0  and  the  surface  temperature  of  the  ice 
(Tg)  is  equal  to  T^,  Eq.  (4)  yields  the  classic  result  that  h  is  a 
simple  function  of  0^/^.  Taking  into  account  the  difference  between  T^ 
and  Tjj  gives  rise  to  the  linear  term  in  (4).  With  no  snow,  the  linear 
term  rapidly  decreases  in  importance  once  h  exceeds  about  20  cm.  A 
small  amount  of  snow,  however,  greatly  extends  the  range  over  which  the 
linear  term  dominates  the  growth. 

With  nominal  values  for  the  thermal  properties  and  h^  -  0,  (4)  can 
be  written  as  h^  +  (13.1hg  +  16.8)h  -  12.96.  Comparison  with  (3)  shows 
that  the  theoretical  coefficients  multiplying  the  h  and  6  terms  are 
substantially  larger  than  the  ones  found  empirically.  If  hg  -  0,  (4) 
predicts  much  too  rapid  a  growth.  With  a  few  centimeters  of  snow  the 
theoretical  growth  is  too  slow  initially  and  too  rapid  once  the  ice 
becomes  thicker.  One  factor  contributing  to  this  problem  is  that, 
contrary  to  what  was  assumed  in  the  derivation  of  (4),  heat  conduction 
in  the  ice  is  not  usually  the  dominant  heat  transfer  mechanism  in  the 
early  stages  of  ice  formation.  Strictly,  then,  (4)  should  be  applied 
only  after  a  solid  ice  cover  has  been  established.  Another  complication 
is  the  fact  that  hg  cannot  be  expected  to  remain  constant  over  extended 
periods  of  time.  It  is,  however,  possible  to  take  into  account  hg(t)  by 
applying  (4)  in  a  stepwise  fashion  between  periods  of  snow  accumulation. 
Such  an  approach  was  taken  by  Nakawo  and  Sinha  (1981)  in  analyzing  ice 
growth  data  from  the  Canadian  Arctic.  Surprisingly,  they  found  that 
growth  predictions  made  by  Including  observed  variations  in  hg  were  less 
accurate  than  those  obtained  using  a  large  constant  snow  depth.  The 
reason  for  the  apparent  failure  of  the  more  rigorous  method  was  probably 
due  to  the  oceanic  heat  flux  at  the  bottom  of  the  ice  (F^) .  Allison 
(1981)  reports  values  of  as  high  as  50  W  m  under  young,  rapidly 
growing  ice,  which  would  be  sufficient  to  decrease  growth  rates  by  1-2 
cm/day . 


Unfortunately,  a  simple  analytic  solution  to  the  heat  balance 
equations  cannot  be  obtained  when  Is  nonzero,  making  It  necessary  to 
employ  numerical  techniques  to  estimate  the  Ice  growth.  In  the  case  of 
young  Ice  the  numerical  procedures  are  relatively  simple  and  growth 
calculations  straightforward.  Using  a  numerical  model  to  describe  the 
thermodynamics  of  young  Ice  (e.g.,  Maykut,  1978)  has  the  advantage  that 
many  of  the  simplifying  assumptions  (e.g.,  zero  radiation  balance)  that 
were  needed  to  obtain  an  analytical  expression  are  unnecessary.  More¬ 
over,  such  a  model  Is  general  In  that  It  Is  not  restricted  to  a 
particular  location  or  season,  provided  that  sufficient  Information  on 
the  thermal  forcing  Is  known.  If  actual  data  on  snowfall,  T^,  F^,  wind 
speed,  and  the  radiative  fluxes  are  available,  very  accurate  growth 
calculations  can  be  made.  Even  If  direct  observations  are  lacking.  It 
Is  possible  to  Improve  on  empirical  predictions  through  the  use  of 
climatological  data  and  suitable  parameterlzatlons  for  the  radiative 
fluxes  {Maykut,  In  press).  The  most  serious  source  of  error  In  such  a 
calculation  would  probably  be  uncertainties  In  F„.  If  only  a  crude 
estimate  of  Ice  growth  Is  needed,  then  the  formula  given  by  Zubov 
(1945),  h^  +  50h  -  80,  will  probably  yield  acceptable  results  as  It 
represents  an  average  of  many  years  of  observation  from  the  Soviet 
Arctic.  Because  It  contains  Implicit  assumptions  regarding  F^(t)  and 
hg(t).  Its  applicability  to  the  Southern  Ocean  Is  doubtful. 

Although  estimating  the  thickness  of  the  Ice  Is  of  primary 
Importance  In  many  practical  problems,  the  rate  at  which  It  grows 
affects  the  level  of  convective  activity  and  the  structure  of  the 
underlying  ocean.  This  Is  because  freezing  removes  freshwater  from  the 
surface  of  the  water  column,  producing  In  effect  a  salt  flux  at  the 
bottom  of  the  Ice.  The  rate  of  Ice  production  Is  extremely  sensitive  to 
Ice  thickness,  dropping  by  roughly  an  order  of  magnitude  In  the  first  50 
cm  of  growth.  Figure  5,  based  on  Anderson's  observations,  shows  growth 
rate  f  as  a  function  of  h  for  several  different  air  temperatures.  Note 
that  the  sensitivity  of  f  to  T^  tends  to  decrease  as  h  Increases. 
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Fig.  5.  Dependence  of  growth  rates  in  young  sea  Ice  on  Ice 

thickness  and  air  temperature  (°C)  (after  Haykut,  in 
press) . 

Because  of  its  low  thermal  conductivity,  snow  Insulates  the  ice  and 
decreases  the  rate  at  which  it  can  lose  heat  to  the  atmosphere.  With  a 
nvimerlcal  model  it  is  possible  to  quantify  the  response  of  young  ice  to 
changes  in  h^  (Fig.  6).  Increasing  h^  decreases  f,  the  effect  being 
much  more  pronoxuiced  for  the  thinner  ice.  For  the  case  illustrated,  f 
is  essentially  constant  if  h^  >  15  cm.  The  reason  for  this  behavior  is 
best  explained  in  terms  of  the  thermal  conductance  (y)  of  the  combined 
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Fig.  6.  Growth  rates  as  a  function  of  snow  depth  for  ice 
thicknesses  of  10  and  50  cm. 

ice -snow  slab.  The  magnitude  of  f  depends  on  the  rate  at  which  heat  is 
conducted  through  the  slab  (F^,) ,  which  for  yoiing  ice  can  be  written 

Fc  -  Y(Tg  -  Tf)  ,  (5) 

where  Tg  is  the  surface  temperature  of  the  slab  and  Y  =  k^kg/Ck^hg  -i- 
kshf) .  Values  of  Y(h,hg)  are  presented  in  Section  4  (see  Fig.  14). 

When  hg  is  small,  y  (and  hence  f)  decreases  rapidly  with  increasing  h; 
when  hg  is  large,  Y  is  small  and  changes  slowly  with  both  h  and  hg. 
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The  principal  simplification  used  in  treating  the  growth  of  young 
ice  was  the  assumption  of  a  linear  temperature  gradient  in  the  ice, 
causing  f  to  respond  immediately  to  changes  in  thermal  forcing  at  the 
surface.  However,  measurements  of  the  temperature  field  in  multiyear 
ice  (Untersteiner,  1961)  indicate  significant  departures  from  linearity, 
even  during  midwinter  when  temperature  fluctuations  are  relatively 
small.  Figure  7  shows  theoretical  predictions  of  f(h)  as  a  function  of 
season  for  thicker  ice  in  the  Central  Arctic.  The  curves  demonstrate 
that  differences  in  thermal  mass  strongly  affect  how  the  ice  responds  to 
conditions  at  the  upper  surface.  In  November,  over  two  months  after 
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Growth  rates  in  multiyear  and  thick  first-year  sea  ice  as 
a  function  of  season  in  the  Central  Arctic  (after  Haykut, 
in  press) . 


freezeup,  ice  thicker  than  about  3  m  continues  to  ablate  because  the 
fall  cooling  has  had  Insufficient  time  to  reach  the  lower  part  of  the 
Ice.  Spring  warming  Is  reflected  In  decreasing  values  of  f  In  the 
thinner  Ice  during  April  and  May,  while  Ice  thicker  than  4  m  continues 
to  feel  the  effects  of  the  winter  cooling.  In  June,  accretion  can  occur 
at  the  bottom  of  the  ice  while  vigorous  melting  Is  taking  place  at  the 
surface.  Clearly  the  assumption  of  a  linear  temperature  profile  Is 
Invalid  for  thicker  Ice  whose  growth  depends  more  on  Its  thermal  history 
than  on  the  Immediate  heat  balance  at  the  surface. 

Heat  transport  and  temperature  changes  In  thicker  ice  can  be 
predicted  by  utilizing  the  standard  second-order  partial  differential 
equation  describing  heat  conduction  In  a  solid.  Combining  this  equation 
with  heat  balance  equations  at  the  upper  and  lower  boundaries  of  the  Ice 
results  In  a  model  that  can  simulate  the  growth  of  any  thickness  of 
undeformed  Ice  (see  Section  SB) .  A  detailed  model  of  this  type  which 
takes  Into  account  effects  of  salinity  and  Internal  solar  heating  has 
been  formulated  and  successfully  tested  under  arctic  conditions  by 
Maykut  and  Vntersteiner  (1971).  A  serious  drawback  to  Incorporating  the 
heat  conduction  equation  Into  an  Ice  model  Is  that  an  accurate  solution 
Involves  strong  numerical  constraints  and  substantial  computer 
resources.  The  problem  can,  to  a  large  extent,  be  overcome  by  only 
crudely  resolving  temperatures  within  the  ice.  Numerical  experiments 
(Semtner,  1976;  Maykut,  In  press)  Indicate  that  results  of  sufficient 
accuracy  for  climate  studies  or  biological  applications  can  be  obtained 
with  as  few  as  three  or  four  Internal  temperature  points.  Numerical 
difficulties  in  this  case  are  minimal. 

E.  Eoulllbrlum  Growth 

The  concept  of  an  equilibrium  thickness  (h^)  Is  frequently  useful 
In  discussing  perennial  Ice  In  the  Arctic.  Such  Ice  loses  mass  each 
summer  by  ablation  at  the  upper  and  lower  boundaries  and  gains  mass  the 
rest  of  the  year  by  accretion  at  the  bottom.  The  amount  of  Ice  accreted 
each  year  depends  strongly  on  h.  If  Interannual  variations  In  atmos¬ 
pheric  and  oceanic  forcing  are  small,  the  ice  will  ultimately  attain  a 
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steady- state  thickness  pattern  where  sumner  ablation  essentially 
balances  the  net  annual  accretion.  The  Ice  Is  then  In  thermodynamic 
equilibrium  with  Its  environment,  and  Ice  thicker  than  equilibrium  will 
experience  net  annual  thinning  while  Ice  thinner  than  equilibrium  will 
undergo  net  growth.  Model  simulations  (Section  5B)  Indicate  that  the 
magnitude  of  h^  Is  especially  sensitive  to  factors  that  affect  the 
amount  of  surface  melting:  air  temperature,  Ice  albedo.  Incoming  radi¬ 
ation,  and  turbulent  heat  exchange.  Despite  small  year  to  year  changes 
In  the  present  climate  and  advectlon  of  Ice  between  different  climatic 
regimes,  observations  show  that  undeformed  older  Ice  In  the  Arctic  Basin 
almost  Invariably  has  a  thickness  In  the  3-4  m  range,  In  agreement  with 
theoretical  expectations  for  h^  In  the  region  (Maykut  and  Vntersteiner , 
1971) .  Althou^  equilibrium  Ice  In  the  Arctic  may  conceivably  exist  for 
any  length  of  time.  It  Is  composed  of  Ice  whose  age  does  not  generally 


exceed  7-10  years.  Over  an  annual  cycle  some  40-50  cm  of  Ice  are  lost 
from  the  surface  and  replaced  by  an  equal  amount  of  new  Ice  at  the 
bottom.  Thus  an  Ice  particle  will  necessarily  move  upward  with  time  at 
a  velocity  of  40-50  cm  per  year  until  eventually  melting  away  at  the 
surface . 

F.  Underwater  Frazil  Production 

To  this  point  the  discussions  have  focused  on  congelation  ice 
(sometimes  called  columnar  Ice)  which  grows  at  the  Ice-water  Interface 
In  response  to  conductive  heat  losses  along  the  temperature  gradient  In 
the  Ice.  It  has  been  widely  accepted  that  such  Ice  makes  up  the  vast 
majority  of  the  polar  Ice  pack;  however,  recent  data  suggest  that  there 
are  some  regions  where  other  Ice  formation  mechanisms  may  be  Important. 
As  mentioned  earlier,  multiyear  Ice  In  the  Weddell  Sea  appears  to  be 
substantially  thicker  than  would  be  expected  on  the  basis  of  congelation 
growth  alone.  Ice  core  data  from  the  region  Indicate  that  congelation 
Ice  accounts  for  only  about  20%  of  the  multiyear  Ice  and  47%  of  the 
first-year  Ice,  with  fine  grained  layers  of  frazil  crystals  making  up 
most  of  the  remainder.  This  means  that  3-4  m  of  frazil  crystals  must 
accumulate  and  be  Incorporated  Into  the  Ice  cover  over  one  or  two  annual 


cycles.  The  higher  percentage  of  frazil  in  thicker  Ice  presumably 
reflects  a  decrease  In  the  rate  of  congelation  growth  relative  to  frazil 
Ice  growth,  suggesting  that  frazil  production  Is  not  strongly  dependent 
on  Ice  thickness.  The  fact  that  frazil  layers  are  often  sandwiched 
between  layers  of  congelation  Ice  Indicates  that  frazil  production  Is 
episodic. 

Underwater  frazil  Ice  also  appears  to  be  common  along  the  antarctic 
coast.  For  example,  Dayton  et  al.  (1969)  observed  1-4  m  thick  billows 
of  frazil  crystals  beneath  the  Ice  In  McMurdo  Sound.  On  some  occasions 
Ice  crystals  appeared  throughout  the  water  column  accompanied  by  the 
formation  of  anchor  Ice  on  the  sea  bottom  and  on  lines  hanging  In  the 
water.  Similar  conditions  were  also  noted  by  Morecki  (1965)  near  Mirny. 
It  Is  suspected  that  periodic  Incursions  of  supercooled  water  generated 
at  nearby  Ice  shelves  may  be  responsible  for  the  coastal  observations, 
but  It  does  not  seem  likely  that  this  mechanism  could  explain  the 
extensive  distribution  of  frazil  In  the  Weddell  Sea. 

Because  frazil  crystals  usually  nucleate  at  some  depth  within  the 
water  column  and  then  float  to  the  surface,  frazil  accumulations  tend  to 
concentrate  sediment  and  algae  contained  in  the  water.  Weeks  and  Ackley 
(1982)  mention  that  concentrations  of  algae  In  young  Weddell  Sea  Ice  can 
be  several  times  that  In  the  water.  Microscopic  examination  of  the  Ice 
revealed  considerable  numbers  of  cells  within  Individual  crystals, 
suggesting  that  such  cells  are  an  effective  nucleating  agent  for  frazil 
Ice.  At  the  same  time,  one  to  two  orders  of  magnitude  more  cells  were 
located  at  grain  boundaries  than  within  crystals.  Since  the  ice  was  too 
young  for  significant  biological  growth  to  have  taken  place  within  the 
Ice,  It  follows  that  the  algae  loading  was  primarily  the  result  of 
scavenging  by  the  rising  crystals. 

Analogous  phenomena  also  appear  to  be  widespread  In  coastal 
portions  of  the  Arctic  where  divers  frequently  encounter  accumulations 
of  slush  on  the  underside  of  congelation  ice.  This  "slush  Ice"  can 
reach  thicknesses  of  at  least  2.5  m,  but  over  50%  of  the  layer  appears 
to  consist  of  unfrozen  water  and  divers  can  easily  reach  Into  It  with 
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their  hands  (Larsen,  1980).  Although  frazil  crystals  form  the  bulk  of 
this  slush,  large  numbers  of  silt  particles  are  often  located  In 
Interstices  between  the  frazil  crystals.  The  source  of  many  of  these 
particles  Is  likely  to  be  scavenging  from  the  water  column,  but  direct 
entrainment  by  Ice  crystals  forming  on  the  sea  floor  could  also  be  a 
factor.  With  time,  upward  loss  of  heat  through  the  overlying 
congelation  Ice  will  cause  Interstitial  water  In  the  slush  to  freeze, 
producing  a  solid  layer.  Normal  congelation  Ice  may  then  form  below  the 
frozen  slush.  This  was  clearly  the  case  In  many  of  the  Weddell  Sea 
cores.  Biological  activity  Is  severely  Impacted  by  the  presence  of 
silt-laden  ice,  light  transmission  and  photosynthesis  being  negligible 
beneath  It.  Moreover,  advectlon  and  subsequent  melting  of  such  Ice 
after  the  spring  breakup  could  represent  an  Important  sediment  transport 
mechanism  In  the  Arctic.  Because  most  of  the  underlce  diving  has  taken 
place  within  a  few  kilometers  of  the  coast.  It  Is  not  known  whether  such 
Intense  frazil  generation  also  occurs  In  deeper  water  where  little  If 
any  sediment  would  be  entrained.  In  the  Central  Arctic  frazil 
production  takes  place  In  leads  and  In  water  below  the  Ice,  but  there  is 
little  evidence  from  core  data  to  suggest  a  significant  effect  on  Ice 
thickness . 

Numerous  mechanisms  have  been  proposed  to  explain  the  occurrence  of 
frazil  Ice  In  the  ocean  (Larsen,  1980;  Weeks  and  Ackley,  1982;  Sewbury, 
1983) .  All  depend  on  having  sllgjit  supercooling  and  some  degree  of 
turbulence  In  the  water  column.  Wind  Induced  mixing  Is  one  such 
mechanism.  However,  Its  effects  are  only  felt  near  the  surface  and  are 
rapidly  damped  out  as  the  grease  ice  accumulates.  In  large  leads  and 
polynyas,  wind  and  wave  action  herd  the  frazil  crystals  to  the  downwind 
side  where  It  can  accumulate  to  depths  of  up  to  1  m  (Martin,  1981). 
Currents  and  local  circulations  can  carry  some  of  these  crystals  beneath 
the  surrounding  Ice  but,  except  for  regions  with  semipermanent  polynyas, 
the  process  Is  self-limiting  and  does  not  seem  capable  of  explaining  the 
widespread  distribution  of  frazil  in  areas  like  the  Weddell  Sea. 

Another  source  of  underwater  frazil  Ice  is  the  upward  movement  of 
water  at  Its  freezing  point.  Decreasing  pressure  as  the  water  rises 
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raises  the  freezing  point,  producing  supercooling  and  the  potential  for 
frazil  formation.  This  mechanism  Is  believed  to  be  responsible  for 
frazil  generation  near  Ice  shelves.  Icebergs,  and  perhaps  even  pressure 
ridge  keels.  Likewise,  contact  between  water  masses  of  different 
salinity,  but  which  are  each  at  their  freezing  point,  can  cause  Ice  to 
form  at  their  Interface  through  the  process  of  double  diffusion.  Such 
Ice  Is  known  to  occur  where  rivers  empty  Into  the  sea  and  during  summer 
runoff  of  meltwater  from  the  Ice  pack.  In  some  areas  like  the  Beaufort 
Sea  coast,  permafrost  extends  beneath  the  sea  floor,  promoting  downward 
conduction  of  heat  and  the  formation  of  anchor  Ice.  Periodic  release  of 
anchor  Ice  could  thus  contribute  to  sediment  loading  and  Ice  accumu¬ 
lation  at  the  bottom  of  the  Ice  cover.  Vhlle  most  of  these  mechanisms 
are  restricted  to  very  special  situations,  underwater  frazil  Ice  can 
also  form  as  a  result  of  thermohallne  convection,  a  process  associated 
with  sea  Ice  throughout  the  polar  regions.  Drainage  of  cold  dense  brine 
from  the  Interior  of  sea  Ice  occurs  during  much  of  the  year  and  produces 
descending  plumes  of  brine  that  cause  supercooling  In  the  adjacent 
water.  Some  of  the  resulting  Ice  forms  as  hollow  stalactites  around  the 
brine  where  It  exits  the  Ice,  but  laboratoiry  studies  Indicate  that  as 
much  as  half  of  the  Ice  production  Is  In  the  form  of  frazil  crystals 
{Martin,  1974).  Similar  plumes  also  exist  under  refreezlng  leads. 

Weeks  and  Ackley  (1982)  feel  that  there  Is  a  greater  potential  for 
frazil  production  In  the  Antarctic  than  In  the  Arctic  and  suggest  that 
thermohallne  convection  could  account  for  much  of  the  frazil  observed  In 
the  Weddell  Sea.  None  of  these  processes  have  yet  been  studied  In 
sufficient  detail  to  allow  quantitative  calculations.  It  appears, 
however,  that  underwater  frazil  Ice  Is  much  more  widespread  than 
previously  thought  and  that  no  single  mechanism  Is  likely  to  explain  all 
of  the  observed  phenomena. 

G.  The  Summer  Melt  Cycle 

Existing  data  Indicate  significant  differences  In  summer  melt 
progression  between  regions  of  perennial  Ice,  seasonal  Ice,  and  coastal 
Ice,  and  between  the  Arctic  and  Antarctic.  In  the  Arctic,  melting  of  the 


snow  cover  begins  as  the  air  temperature  approaches  the  freezing  point 
and  typically  takes  2-3  weeks.  Accompanying  the  snow  melt  is  a  dramatic 
decrease  in  albedo  due  to  the  puddling  of  melt  water  on  the  surface . 

Melt  pond  coverage  on  multiyear  ice  may  briefly  reach  values  of  50-60%, 
but  quickly  decreases  to  about  30%  once  the  snow  cover  has  been  removed. 
With  the  development  of  melt  ponds,  area- averaged  albedos  drop  from 
about  0.7  to  0.5.  Deepening  and  consolidation  cause  a  continual 
decrease  in  the  areal  extent  of  the  ponds  throughout  the  summer.  Values 
as  low  as  10%  are  common  by  the  end  of  August.  Albedos  increase  from 
about  0.50  to  0.55  during  this  period.  Because  turbulent  heat  fluxes 
and  longwave  radiation  act  to  remove  heat  from  the  surface  during  the 
summer,  the  energy  needed  to  drive  the  melt  cycle  in  regions  of 
perennial  arctic  ice  comes  almost  entirely  from  shortwave  radiation. 

Initial  melting  on  seasonal  ice  produces  nearly  complete  pond 
coverage  because  of  the  lack  of  surface  topography.  Even  when  the  melt 
cycle  is  well  advanced,  it  is  not  unusual  for  shallow  ponds  to  occupy 
60%  or  more  of  the  surface  and  average  albedos  to  be  in  the  0.40-0.45 
range.  In  coastal  regions  dust  deposited  from  nearby  land  areas  tends 
to  lower  the  albedos  even  further.  Areally  averaged  albedos  as  low  as 
20%  have  been  reported  from  such  regions  {Langleben ,  1966).  Air 
temperatures  in  coastal  areas  are  typically  warmer  than  the  ice,  and  the 
turbulent  fluxes  make  a  substantial  contribution  to  the  melting.  Lower 
albedos  and  thinner  ice  also  allow  greater  input  of  solar  energy  to  the 
underlying  ocean  in  areas  of  seasonal  ice.  Thus,  even  when  the  Incoming 
radiation  fluxes  are  similar,  ice  decay  proceeds  more  rapidly  in  coastal 
regions  and  areas  of  first-year  ice  than  in  the  interior  pack. 

The  pattern  of  ice  decay  in  the  Southern  Ocean  appears  to  differ 
strikingly  from  that  in  the  Arctic.  Katabatic  winds,  higher  latitudes, 
and  lower  oceanic  heat  fluxes  cause  melt  rates  near  the  continental 
boundaries  to  be  small  compared  with  those  near  the  free  boundary  of  the 
antarctic  ice  pack.  Moreover,  the  melt  ponds  which  are  so  character¬ 
istic  of  melting  ice  in  the  Arctic  have  rarely  been  observed  in  the 
Antarctic,  and  surface  melting  appears  to  be  minimal.  The  reason 
appears  to  be  related  to  differences  in  the  meteorological  variables 


that  control  the  surface  heat  balance.  Surface  melting  in  the  Arctic 
often  takes  place  with  air  temperatures  significantly  below  freezing. 

In  the  Antarctic,  however,  Andreas  and  Ackley  (1982)  argue  that  lower 
relative  humidities  and  stronger  winds  enhance  turbulent  heat  losses 
from  the  Ice  and  allow  surface  melting  to  occur  only  with  air  tempera¬ 
tures  significantly  above  0°C,  a  condition  rarely  encountered  over  the 
Ice  cover  In  the  Southern  Ocean.  In  spite  of  the  near  absence  of 
surface  melting,  most  of  the  antarctic  Ice  pack  disappears  each  summer, 
presumably  as  a  result  of  melting  at  the  Ice-water  Interface.  Gordon 
(1981)  estimates  that  no  more  than  half  of  the  necessary  energy  could  be 
derived  from  the  relatively  warm  deep  water  below  the  Southern  Ocean 
pycnocllne.  The  remainder  must  come  from  shortwave  radiation  absorbed 
by  the  upper  ocean,  either  through  leads  In  the  Ice  pack  or  near  the  Ice 
edge. 

An  open  lead  absorbs  over  90%  of  the  Incident  shortwave  radiation, 
about  half  of  which  Is  absorbed  In  the  upper  two  meters  of  the  water. 
Part  of  this  energy  goes  to  lateral  melting  on  floe  edges,  causing  an 
Increase  In  the  area  of  open  water  (A^) .  This  represents  a  positive 
feedback  process  which  affects  both  the  rate  of  retreat  of  the  Ice  edge 
and  the  overall  heat  and  mass  balance  of  the  Ice  pack.  In  sheltered 
bays  and  fjords  along  the  arctic  coast,  for  example,  surface  melting  Is 
generally  Inadequate  to  remove  the  2-2.5  m  of  Ice  which  grows  each 
winter.  The  annual  disappearance  of  this  Ice  can  only  be  explained  by 
taking  Into  account  the  accelerated  decay  produced  by  solar  heating  In 
leads  (Langleben,  1972). 

Both  Ice  dynamics  (diverging  motions)  and  thermodynamics  (lateral 
melting  on  floe  edges)  affect  A^,  though  field  data  are  not  yet  adequate 
to  Identify  their  relative  Importance  In  the  different  regions.  Except 
In  areas  like  the  Greenland  Sea,  Ice  divergence  Is  believed  to  be  much 
smaller  In  the  Arctic  than  In  the  Antarctic.  Nevertheless,  substantial 
amounts  of  open  water  exist  In  the  Arctic  where  the  summer  Ice  cover 
undergoes  a  dramatic  Increase  In  the  number  and  size  of  leads.  Even  In 
regions  of  perennial  Ice,  summer  values  of  A,,  can  reach  10-15%  (Hall, 
1978).  Clearly  the  effects  of  A^(t)  on  the  Input  of  solar  energy  and 


meltwater  to  the  upper  ocean  can  have  a  substantial  Impact  on  biological 
processes  there. 

Various  attempts  have  been  made  to  relate  the  summer  decrease  in 
ice  thickness  to  through  empirical  equations  analogous  to  those 
developed  for  young  growing  Ice.  Bilello  (1961),  for  example,  found  a 
correlation  of  0.93  between  his  data  and  thickness  changes  given  by  the 
equation 

Ah  -  0.550'  ,  (6) 

where  Ah  -  the  total  decrease  in  Ice  thickness  and  6 '  -  the  accumulated 
degree-days  above  freezing.  Bilello  also  mentions  a  similar  relation¬ 
ship  [Ah  -  0.51(0"  -  32),  where  0"  B  accumulated  degree-days  above  -5°C] 
developed  by  Karelin  for  the  Soviet  Arctic.  The  high  degree  of  corre¬ 
lation  between  and  Ah  Is  surprising  because  shortwave  radiation  Is 
the  dominant  term  In  the  summer  energy  balance.  Variations  in  cloud 
cover,  surface  albedo,  ice  concentration,  and  would  seem  to  be  as 
Important  as  T^.  The  probable  explanation  Is  that  the  observations  were 
taken  close  to  shore  and  were  strongly  Influenced  by  warm  air  from 
nearby  land  masses.  It  is  unlikely  that  such  relationships  could  be 
employed  away  from  coastal  areas  where  summer  values  of  T^  are 
invariably  close  to  the  freezing  point. 


3.  STRUCTURE  AND  SALINITY 
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A.  Crystal  Orientation 


Ice  crystals  are  hexagonal  in  form  and  are  characterized  by  one 
principal  hexagonal  axis  of  symmetry  (the  c-axis)  which  is  the  optic 
axis  of  the  crystal.  Perpendicular  to  the  c-axls  is  the  basal  (0001) 
plane  defined  by  three  crystallographic  axes  (the  a-axes)  of  equal 
length  with  angles  of  120°  between  the  positive  ends.  Six  crystal  faces 
parallel  the  c-axis  and  are  often  called  prism  faces.  The  arrangement 
of  atoms  in  an  ice  lattice  is  such  that  fracture  along  the  basal  plane 
requires  the  breaking  of  only  two  bonds,  while  fracture  along  any  plane 
normal  to  this  requires  the  breaking  of  four  or  more  bonds.  Slip  and 
cleavage  are  thus  most  likely  to  occur  parallel  to  the  basal  plane. 

In  natural  sea  ice  the  orientation  of  individual  crystals  tends  to 
be  highly  organized.  Only  near  the  surface  where  the  ice  is  composed  of 
frozen  slush  or  grease  ice  are  crystal  orientations  random.  Once 
congelation  growth  begins,  crystals  whose  c-axis  are  parallel  to  the 
ice -water  interface  quickly  begin  to  dominate  the  structure  of  the  ice 
sheet.  The  explanation  for  this  phenomenon  is  anisotropic  growth  on 
different  crystal  faces  which  results  in  a  geometric  selection  of 
crystals  with  favorable  orientations.  The  structure  of  the  ice  lattice 
is  such  that  the  density  of  potential  bonding  sites  is  greater  on  a 
prism  face  than  on  a  basal  plane.  This  means  that,  for  a  given  degree 
of  supercooling,  growth  parallel  to  the  c-axls  will  be  much  slower  than 
growth  perpendicular  to  the  c-axls.  Observations  (Hallett,  1960)  show 
growth  rate  differences  of  as  much  as  two  orders  of  magnitude.  Crystals 
whose  c-axes  are  closer  to  the  horizontal  grow  downward  into  the  water 
faster  than  those  with  a  more  vertical  orientation,  wedging  out  the 
slower  growing  crystals.  This  process  proceeds  quite  rapidly  until  only 
those  crystals  with  c-axes  parallel  to  the  freezing  Interface  remain. 

As  we  shall  see  later,  this  turns  out  to  be  an  important  factor  in 
determining  the  structure  and  properties  of  sea  ice. 


The  region  over  which  the  crystal  orientation  shifts  from  being 
random  to  having  essentially  all  c-axes  horizontal  is  termed  the 
transition  zone  and  typically  occurs  over  a  vertical  distance  of  only 
5>10  cm.  The  transition  zone  begins  at  the  bottom  of  the  frozen  frazil 
layer  and  hence  may  be  located  anywhere  in  the  upper  meter  of  the  ice 
cover,  depending  on  the  initial  thickness  of  the  frazil  layer.  Beneath 
the  transition  layer  there  is  little  interference  between  the  growth  of 
adjacent  crystals  and  continued  development  of  these  crystals  occurs 
primarily  in  the  vertical.  The  region  between  the  bottom  of  the 
transition  zone  and  the  ice-water  Interface  is  the  so-called  columnar 
zone  and  is  characterized  by  crystals  strongly  elongated  in  the 
direction  of  the  temperature  gradient.  Crystal  diameter  in  the  upper 
half  meter  of  the  columnar  zone  is  typically  about  1-2  cm.  There  is  a 
definite  tendency  for  average  crystal  diameter  to  increase  with  depth  in 
a  roughly  linear  fashion  (Weeks  and  Hamilton,  1962;  Weeks  and  Ackley, 
1982).  The  reason  appears  to  be  associated  with  the  decreasing  growth 
rate  of  the  ice  as  it  thickens. 

Until  recently  it  was  generally  accepted  that  horizontal  c-axis 
alignments  were  random  throughout  the  columnar  zone.  However,  several 
studies  (e.g.,  Cherepanov,  1971;  Weeks  and  Cow,  1978)  now  indicate  that 
strong  horizontal  alignments  exist  over  large  portions  of  the  Arctic. 


Near  shore  observations  suggest  that  the  c-axes  are  closely  oriented  in 
the  direction  of  the  average  current  at  the  bottom  of  the  ice.  Since 
the  microscale  roughness  on  the  underside  of  the  ice  is  greater  when  the 
current  is  parallel  to  the  c-axis  (see  next  section) ,  Weeks  and  Cow 
(1978)  argue  that  the  resulting  increase  in  turbulence  decreases  the 
thickness  of  the  stable  boundary  layer  at  the  growing  interface,  giving 
such  crystals  a  slight  growth  advantage.  Although  this  explanation  is 
plausible,  additional  field  and  laboratory  work  is  needed  to  verify  it. 

B.  The  Skeletal  Laver  and  Brine  Entrapment 

Examination  of  the  underside  of  growing  sea  ice  reveals  the 
presence  of  a  delicate  layer  of  thin  ice  platelets  (Fig.  8a)  arranged 
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Fig.  8a.  Schematic  illustration  of  the  platelet  structure  on  the 
bottom  of  a  growing  ice  crystal. 


perpendicular  to  the  c-axes  of  the  individual  crystals.  The  many 
platelets  growing  on  the  bottom  of  each  exposed  crystal  are  generally 
parallel  (Fig.  8b)  with  a  spacing  (a^)  of  0. 5-1.0  mm.  The  zone  of 
platelet  occurrence  has  a  vertical  dimension  of  about  1-3  cm  and  is 
frequently  referred  to  as  the  skeletal  layer.  It  has  been  found  both 
theoretically  and  experimentally  that  a^  is  proportional  to  . 

Clearly,  flow  parallel  to  the  rows  of  platelets  generates  the  least 
turbulence,  allowing  a  slightly  thicker  boundary  layer  to  develop 
beneath  crystals  with  c-axes  normal  to  the  flow  and  possibly  explaining 
the  apparent  growth  advantage  of  crystals  more  closely  aligned  with  the 
flow  (Weeks  and  Gov,  1978). 
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Fig.  8b.  Plan  view  showing  platelet  arrangement  at  the  underside 
of  growing  sea  Ice.  Curved  lines  represent  crystal 
boundaries  and  short  straight  lines  the  platelets.  Note 
that  platelets  on  an  individual  crystal  tend  to  form  In 
parallel  sets  normal  to  the  horizontal  direction  of  the 
c-axls  (double  ended  arrow) . 

Absence  of  a  skeletal  layer  beneath  freshwater  Ice  suggests  that 
the  complex  microstructure  found  at  the  growing  interface  of  sea  Ice  Is 
somehow  related  to  the  presence  of  salt  In  the  water.  When  seawater 
freezes,  most  of  the  impurities  (salts)  are  rejected  from  the  Ice 
lattice,  resulting  in  the  formation  of  a  thin  boundary  layer  of  very 
salty  water  ahead  of  the  advancing  Interface.  The  continual  addition  of 
salt  causes  the  salinity  at  the  top  of  this  layer  to  be  greater  than  the 
salinity  beneath  the  layer.  Because  of  the  higher  salinity,  the 
temperature  at  the  Ice-water  Interface  (T^)  will  be  lower  than  that  In 
the  underlying  water.  These  conditions  produce  a  downward  flux  of  salt 
and  an  upward  flux  of  heat  through  the  layer.  If  we  assume  that  the 
transport  of  heat  and  salt  Is  primarily  the  result  of  diffusive 
processes.  It  follows  that  constitutional  supercooling  can  occur  In  this 
boundary  layer.  Because  of  large  differences  In  the  rate  at  which  heat 
and  salt  diffuse,  the  flux  of  heat  through  the  boundary  layer  Is  much 
more  rapid  than  the  flux  of  salt,  so  that  salinity  undergoes  an 
exponential  type  decrease  below  the  Interface  whereas  the  temperature 
Increase  Is  more  nearly  linear.  Since  it  is  the  salinity  that 
determines  the  freezing  point  of  the  solution,  the  actual  temperature 
will  typically  be  below  the  freezing  point  (Fig.  9).  Under  such 
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Fig.  9.  Constitutional  supercooling  in  the  boundary  layer  ahead 
of  an  advancing  ice  Interface. 

conditions,  laboratory  observations  {Harrison  and  Tiller,  1963)  reveal 
that  the  ice -water  Interface  Initially  develops  an  array  of  parallel 
knife-edged  cells  on  which  scallops  form  and  deepen  until  the  cells 
separate  into  a  line  of  individual  platelets.  The  length  and  spacing  of 
these  platelets  act  to  minimize  the  amount  of  supercooling  ahead  of  the 
Interface.  An  approximate  theory  has  been  worked  out  to  describe 
criteria  for  constitutional  supercooling  and  the  geometry  of  the 
Interface  (Weeks  and  Ackley,  1982).  Efforts  are  under  way  to  develop  a 
more  exact  theory  by  allowing  for  free  convection  near  the  interface 
Instead  of  only  diffusion. 

The  significance  of  the  skeletal  layer  is  that  it  provides  a 
mechanism  for  including  substantial  amounts  of  salt  in  the  ice.  As  the 
platelets  lengthen,  they  grow  wider  and  periodically  small  ice  bridges 
form  between  adjacent  platelets,  entrapping  pockets  of  brine.  These 
brine  pockets  are  typically  long  and  narrow,  being  on  the  order  of 
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0.05  mm  In  diameter.  Sea  Ice  structure  Is  thus  characterized  not  only 
by  the  columnar  crystals,  but  also  by  many  vertical  strings  of  liquid 
Inclusions  that  mark  the  original  boundaries  between  the  platelets.  As 
discussed  In  Section  4,  these  brine  pockets  dictate  the  bulk  thermal  and 
mechanical  properties  of  the  Ice.  Without  this  entrapment  mechanism, 
there  would  be  little  to  distinguish  sea  Ice  from  freshwater  Ice. 

Initially,  the  composition  of  the  entrapped  brine  should  have  the 
same  relative  concentration  of  Ions  as  normal  seawater.  However,  as 
temperature  decreases  In  the  surrounding  Ice,  solid  salts  begin  to 
precipitate  out  of  the  brine  solution.  While  mlrablllte  (Na2SO^'10H20) 
begins  to  form  at  -8.2°C,  the  most  Important  precipitate  Is  sodium 
chloride  (NaCil‘2H20)  which  occurs  below  a  temperature  of  -22.9°C. 
Appearance  of  the  NaC£‘2H20  causes  the  Ice  to  take  on  a  milky  color  and 
alters  Its  optical  properties.  Other  solid  salts  that  may  be  present  In 
small  amounts  Include  CaC03‘6H20,  MgC£.2‘8H20,  KCl ,  HgCl2 '  I-2H2O,  and 
CaCJl2'^H2^  (Weeks  and  Ackley,  1982). 

In  general,  the  salinity  of  the  brine  In  a  pocket  Is  Independent  of 
the  salinity  of  the  original  water  because  of  mass  changes  taking  place 
at  the  walls  of  the  pocket.  The  temperature  of  the  brine  (T^) ,  which  In 
the  normal  state  of  thermal  equilibrium  Is  the  same  as  that  of  the 
surrounding  Ice,  must  be  at  the  freezing  point  of  the  solution  (T£) .  If 
>  Tf,  melting  will  occur  on  the  Ice  walls,  lowering  the  salinity  of 
the  brine  (S^)  until  T£  -  Similarly,  If  <  T£,  freezing  will 

occur  In  the  brine,  raising  Sj,  until  T£  and  coincide.  This  being  the 
case,  salinity  In  a  brine  pocket  can  be  related  to  T^  through  the 
following  freezing  point  equations  which  are  based  on  the  work  of  Assur 
(1958): 


1000  Tb 


~  Tb  -  54.11  ’ 


Tb  >  -8.2®C 


(7) 


and 


62.4  -  10.31  Tb 


1.0624  -  0.01031  Tv 


-22.9  <  Tb  4  -8.2®C 


(8) 
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where  is  In  ^/OO  and  in  °C.  Also  fitting  the  data  reasonably  well 
are  the  approximate  relations  given  by  Morey  et  el.  (1984):  -  9.65 

-  14.8  when  -8.2  ^  i  -2°C,  and  -  78.11  -  6.6  Tj,  when  -22.9  i 
^  -8.2°C.  The  change  in  the  relationship  below  -8°C  is  due  to 

the  precipitation  of  Na2SO^'10H2O  from  the  solution.  Below  -23°C  the 
data  are  somewhat  ambiguous,  and  proper  freezing  point  relations  have 
not  definitely  been  established.  Since  the  temperature  of  the 
surrounding  ice  (T^)  will  normally  equal  T^j,  it  is  possible  to 
calculate  by  simply  measuring  and  applying  (7)  or  (8). 

C.  Ice  Salinity 

The  salt  content  of  sea  ice  is  usually  described  in  terms  of  a  bulk 
salinity  (S^)  defined  as 

S  =  -  mass  of  salt_  —  .  ^^3 

^  mass  of  ice  +  mass  of  brine 

where  Sj^  is  in  °/oo.  The  amount  of  salt  initially  entrapped  by  the 
growing  ice  depends  strongly  on  growth  rate  and  water  salinity. 
Laboratory  measurements  by  Cox  and  Weeks  (1975)  show  two  different 
regimes.  When  f  <  2’10*^cm  s'^  ('''1.7  cm  d*^), 

Sj^  -  S^(0.8439  +  0.0529  In  f )  ,  (10) 

while  for  f  >  2 ’10'^  cm  s*^ 


S^  -  0.26  S^[0.26  +  0.74  exp(-7243f)] 


The  flux  of  salt  to  the  water  beneath  the  growing  ice  can  thus  be 
written  as  ”  ^l^^^w  *  ^i^  estimated  using  (10)  or  (11). 

If  the  brine  did  not  migrate  in  the  ice,  we  would  expect  to  see 
stationary  salinity  profiles  with  a  maxlmvim  near  the  surface  and  a 
continuous  decrease  with  depth,  reflecting  the  reduction  in  f  as  the  ice 
thickened.  Instead  we  find  that  S£(z)  is  highly  time- dependent  with  a 


vertical  scatter,  presumably  correlated  with  changes  in  thermal  forcing 
(see  Nakawo  and  Sinha,  1981).  Likewise,  the  magnitude  of  can  be 
quite  different  for  Ice  of  equal  thickness  because  of  differences  In 
growth  rate.  Nevertheless,  Fig.  10  Illustrates  several  Important 
points.  Profiles  In  thinner  Ice  tend  to  be  C-shaped  Indicating  salt 
enrichment  In  the  lower  half  of  the  slab.  The  average  salinity  also 
decreases  as  the  Ice  ages  (l.e.,  becomes  thicker),  suggesting  loss  of 
brine  to  the  ocean.  Finally,  there  Is  a  dramatic  change  In  Sj^(z)  In 
multiyear  and  melting  first-year  Ice. 

To  explain  the  observed  profiles,  Unterstelner  (1968)  examined  the 
role  of  four  mechanisms  that  transport  brine  within  the  Ice.  Brine 
pocket  migration  will  occur  In  the  presence  of  a  temperature  gradient 
across  the  slab.  Melting  on  the  warm  side  of  the  pocket  and  Ice 
accretion  on  the  cool  side  cause  the  pocket  to  move  opposite  to  the 
direction  of  heat  flow.  Under  normal  conditions,  this  mechanism  Is  too 
slow  to  account  for  a  significant  amount  of  brine  transport,  although  If 
the  pocket  Is  large  enough  to  support  convective  heat  transfer,  movement 
can  be  quite  rapid.  Another  mechanism  Is  brine  expulsion  resulting  from 
the  pressure  buildup  In  refreezlng  pockets.  Again  brine  expulsion  does 
not  appear  to  play  a  major  role  In  the  desalination  of  sea  Ice.  It  Is 
likely  to  be  Important  only  during  the  early  stages  of  growth  when  the 
Ice  Is  undergoing  rapid  cooling.  The  "salt  flowers"  and  extremely  high 
salinities  often  found  at  the  surface  of  young  Ice  are  probably  the 
result  of  this  process. 

The  mechanism  primarily  responsible  for  observed  salinity  changes 
In  first-year  Ice  Is  almost  certainly  gravity  drainage.  Because  brine 
pockets  typically  form  In  vertical  strings,  they  tend  to  become 
Interconnected  In  response  to  temperature  changes  and  Internal  stresses, 
forming  long  vertical  tubes  or  brine  channels  which  allow  the  vertical 
movement  of  brine  throu^  the  Ice.  During  most  of  the  year  temperature 
Increases  with  depth,  and  colder  denser  brine  In  the  upper  part  of  the 
Ice  will  sink  downward  to  displace  the  warmer  less  dense  brine.  In 
addition,  bottom  growth  causes  the  Ice  sheet  to  move  upward,  creating  a 
hydrostatic  head  that  forces  brine  downward  In  the  channels  and  out  of 


the  Ice.  Average  brine  channel  diameters  near  the  bottom  of  the  ice  are 
on  the  order  of  0.4  cm.  Depending  on  ice  thickness,  as  many  as  50-300 
channels  per  square  meter  may  exist  along  the  bottom  of  the  ice.  The 
cold  streamers  of  brine  leaving  the  ice  cause  hollow  tubes  of  ice  to 
form  around  them.  These  "ice  stalactites"  have  been  observed  to  reach 
lengths  of  6  m  in  the  Antarctic,  but  lengths  of  25-50  cm  seem  to  be  more 
representative  in  the  Arctic.  Brine  movement  within  the  ice  causes 
complex  changes  within  the  channels.  It  has  been  noted  (Lake  and  Lewis, 
1970)  that  primary  brine  channels  are  often  accompanied  by  smaller 
tributary  channels  which  provide  additional  brine  and  form  what  is,  in 
effect,  a  brine  drainage  system.  Stalactite  formation  and  flow  in  the 
primary  channels  have  been  extensively  studied  (Martin,  1974;  Eide  and 
Martin,  1975;  Niedrauer  and  Martin,  1979),  but  few  details  are  available 
regarding  the  formation  and  role  of  the  tributary  channels. 

The  final  mechanism  is  flushing  by  surface  meltwater.  As  surface 
temperatures  approach  the  freezing  point  in  spring  and  early  summer, 
brine  volume  Increases  and  the  ice  becomes  extremely  porous.  Melting  of 
the  snow  cover  and  the  formation  of  ponds  on  the  surface  allow 
freshwater  to  percolate  into  the  ice,  resulting  in  a  further  decrease  in 
S^.  Calculations  ((Inters teiner,  1968)  indicate  that  the  greatest  salt 
loss  occurs  where  brine  volume  is  a  maximum  and  that  flushing  appears  to 
be  sufficiently  vigorous  in  the  Arctic  to  explain  the  observed 
differences  between  salinity  profiles  in  first-year  and  multiyear  ice. 
Flushing  does  not  appear  to  be  a  significant  factor  in  the  evolution  of 
antarctic  sea  ice. 

Data  from  a  wide  variety  of  sources  and  locations  were  compiled  by 
Cox  and  Weeks  (1974)  to  produce  plots  of  average  salinity  versus  ice 
thickness  in  the  Arctic  (Fig.  11).  For  cold  ice  the  plots  show  a 
roughly  linear  decrease  in  up  to  a  thickness  of  about  40  cm,  at  which 

point  the  slope  abruptly  changes,  possibly  indicating  a  shift  in  the 
dominant  desalination  mechanism  from  brine  expulsion  to  gravity 
drainage.  It  should  be  noted  that,  in  spite  of  a  wide  range  of  growth 
conditions,  ^^(h)  shows  surprisingly  little  scatter.  Measurements  of 


when  the  ice  is  melting  are  complicated  by  brine  leakage  from  the  cores 
but  there  Is  at  best  only  a  weak  dependence  on  h. 

Chemical  studies  (Reeburgh  and  Springer- Young ,  1983)  show  that 
sulfate-chlorlnlty  ratios  In  sea  Ice  are  different  from  the  seawater 
ratio.  On  aging,  sulfate  In  the  Ice  Is  mobilized  and  removed  In  a 
constant  ratio  to  chlorlnlty.  Indicating  conservative  behavior. 
Comparison  of  first-year  and  multiyear  Ice  suggests  that  fractionation 
occurs  during  Ice  formation,  but  that  further  fractionation  within  the 
Ice  does  not  appear  to  take  place.  Because  much  of  the  entrapped  brine 
ultimately  returns  to  the  ocean,  there  seems  to  be  no  lasting  effect  on 
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Fig.  11.  Average  salinity  of  arctic  sea  Ice  as  a  function  of  Ice 
^hlckness  (after  Cox  and  Veeks,  1974).  Soll^  lines  show 
S^(h)  for  cold  Ice,  while  dashed  line  shows  S^Ch)  for 
warm  Ice  late  In  the  melt  season.  Data  points  for  the 
cold  Ice  are  given  by  dots  and  for  the  warm  Ice  by  open 
circles . 


the  chemical  composition  of  the  underlying  water  as  a  result  of  the 
freezing  process. 

The  evolution  of  Ice  salinity  can  be  summarized  as  follows.  High 
Initial  growth  rates  entrap  substantial  amounts  of  brine,  producing 
values  near  the  surface  of  10-20  °/oo.  As  h  Increases,  f  decreases  and 
less  brine  is  entrapped.  However,  cooling  of  the  upper  part  of  the  ice 
results  In  brine  expulsion  that  allows  In  the  lower  part  of  the  Ice 
to  remain  hl^.  Substantial  quantities  of  brine  are  lost  to  the 
underlying  ocean  during  this  period.  By  the  time  h  approaches  40  cm, 
brine  channels  appear  to  be  sufficiently  developed  for  gravity  drainage 
to  control  the  rate  of  desalination.  The  Ice  continues  to  lose  salt 
through  this  process  until  the  onset  of  snow  melting  when  flushing  by 
freshwater  occurs.  Beneath  elevated,  drained  areas  there  1^  a  sharp 
drop  In  salinity  whereas  beneath  low  areas  close  to  the  freeboard  level 
the  decrease  Is  much  less  pronounced.  This  pattern  persists  even  In 
older  multiyear  Ice  where  horizontal  salinity  gradients  exist  which  are 
correlated  with  surface  topography  (Cox  and  Weeks,  1974). 

D.  Annual  Layering  In  Sea  Ice 

Cores  taken  In  multiyear  arctic  sea  Ice  reveal  the  existence  of 
horizontal  layers  which  are  suspected  to  develop  during  the  annual  thaw 
cycles  (Cherepanov,  1957;  Schwarzacher ,  1959).  These  layers  are  of  two 
different  types.  The  first  type  Is  2-5  mm  thick,  milky  white  in  color, 
with  a  sharp  upper  boundary  and  an  Irregular  lower  boundary.  The  origin 
of  this  layer  Is  uncertain,  but  It  appears  to  form  during  the  summer 
after  the  Ice  has  stopped  growing.  It  has  been  suggested  that  the 
mllklness  of  the  Ice  may  be  associated  with  biological  activity  In  the 
ocean  (Weeks  and  Ackley,  1982).  The  second  type  is  much  thicker  (1-10 
cm)  and  appears  to  result  from  the  freezing  of  surface  meltwater  that 
penetrates  beneath  the  Ice.  The  size  of  the  cirystals  Is  much  smaller 
than  In  the  surrounding  layers  of  normal  congelation  Ice  and  the 
salinity  much  lower  (1-1.5  °/oo).  The  presence  beneath  the  ice  of 
freshwater  from  surface  melting  can  be  explained  by  flow  through  natural 


or  artificial  drain  holes  or  by  runoff  into  nearby  leads  that 
subsequently  close.  If  vertical  mixing  is  small,  there  will  be  a  sharp 
interface  between  the  freshwater  at  its  freezing  point  (0®C)  and  the 
underlying  seawater  at  its  freezing  point  ('V'-1.8°C).  This  results  in 
supercooling  in  the  freshwater  near  the  interface  and  ultimately  in  the 
growth  of  ice  crystals  throughout  the  freshwater  layer  (Martin  and 
Kauffman,  1974).  Since  freshwater  does  not  necessarily  occur  beneath 
the  ice  each  summer  at  a  given  location,  the  refrozen  meltwater  layers 
are  often  separated  by  one  or  more  of  the  thin  milky  layers.  These 
layers  make  it  possible  to  estimate  the  age  at  various  depths  within  the 
ice.  Obtaining  the  total  age  of  the  Ice  cover  by  this  method,  however, 
is  difficult  because  melting  and  refreezing  near  the  surface  tend  to 
obliterate  the  uppermost  layers. 


4.  PHYSICAL  PROPERTIES  OF  SEA  ICE 


As  the  ice  changes  temperature,  internal  melting  or  freezing  within 
the  brine  Inclusions  affect  the  brine  volume  of  the  Ice.  The  fractional 
volume  of  brine  In  the  Ice  can  be  written  as 


= 


P^Si 
^  ’ 


(12) 


where  is  the  density  of  the  Ice  and  P^  is  the  density  of  the  brine. 
From  Eqs.  (7)  and  (8)  we  know  that  S^  is  dependent  only  on  Tj^,  and  hence 
will  be  a  function  of  S^  and  T^.  Application  of  (12)  over  a  wide 
temperature  range  Is  complicated  by  the  precipitation  of  solid  salts  and 
uncertainties  In  (^(T^) ;  however,  a  table  of  for  standard  sea  Ice  has 
been  tabulated  by  Assur  (1958)  down  to  a  temperature  of  -54°C.  In  most 
cases  values  of  sufficient  accuracy  can  be  obtained  with  a  simple 
equation  developed  by  Frankenstein  and  Gamer  (1967): 
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vi,  ==  Si(0.0532  -  .  -22.9  <1^5  -0.5'’C  ,  (13) 

where  Sj^  Is  in  °/oo  and  in  percent.  The  brine  volume  at  the  bottom 
of  growing  sea  Ice  can  vary  from  about  8%  to  40%.  For  multiyear  Ice 
Is  typically  about  10%,  except  near  the  surface  during  the  summer  when 
values  as  high  as  40-50%  have  been  observed.  As  will  be  shown  below, 
the  magnitude  of  has  a  major  Impact  on  the  bulk  thermal,  mechanical, 
optical,  and  electromagnetic  properties  of  sea  Ice. 

In  addition  to  the  entrapment  of  brine,  growing  sea  Ice  also 
captures  significant  numbers  of  air  bubbles.  According  to  Schwerdtfeger 
(1963),  the  fractional  volume  of  air  bubbles  In  sea  Ice  (v^)  Is  given  by 
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where  -  0.9178  Mg  m"  is  the  density  of  pure  ice,  is  the  bulk 
density  of  the  sea  ice,  and  -  S^/1000  is  the  fractional  salt  content 
of  the  ice.  Equation  (14)  can  be  rewritten  and  simplified  to  obtain  an 
approximate  relation  for  p^C  v^,T£,Sj^)  , 


0.004563- 

Pi  =  (1  -  Va){l - - - )po  . 


It  can  be  seen  from  (15)  that  p^  is  mainly  determined  by  at  warm 
temperatures  and  by  at  lower  temperatures.  Because  the  major  changes 
in  thermal  and  mechanical  properties  occur  when  the  ice  is  warm,  air 
bubbles  will  be  neglected  in  the  following  discussions. 

A.  Thermal  Properties 

If  we  remove  heat  from  a  block  of  sea  ice,  decreases,  causing 
brine  to  freeze  in  the  pockets.  This  freezing  releases  latent  heat 
which  acts  to  warm  the  ice.  As  a  result,  sea  ice  cools  more  slowly  than 
fresh  ice.  Alternatively,  warming  produces  melting  in  the  brine  pockets 
which  slows  the  rate  at  which  Tj^  increases.  Brine  pockets  thus  behave 
as  thermal  buffers  that  retard  the  transport  of  heat  in  sea  ice,  which 
affects  both  conductivity  and  specific  heat.  These  changes  in  brine 
volume  also  alter  the  latent  heat  of  fusion  of  the  ice. 

1.  Conductivity 

The  conductive  heat  flux  in  the  ice  (F^,)  is  given  by 
Fj,  -  kj^(9T£/3z),  where  k^  is  the  thermal  conductivity  of  the  ice.  The 
magnitude  of  k^  and  the  temperature  gradient  thus  determine  the  rate  of 


conductive  heat  transfer.  In  pure  Ice  the  conductivity  (k^)  depends  on 
temperature  and  is  approximately  -  9.828  exp  (-0.0057T),  where  T  is 
in  °K  and  k^  in  W  m“^  °K'^  (Yen,  1981).  A  typical  value  of  k^  for  warm 
ice  would  be  2.1  W  m'^  (0.005  cal  cm'^  s’^  °K‘^).  There  appears  to 

be  a  slight  difference  between  k^  measured  parallel  to  the  c-axis  and 
that  measured  normal  to  it,  but  this  difference  is  small  (about  5%)  and 
can  usually  be  neglected.  In  sea  ice,  however,  conductivity  is  substan¬ 
tially  lower  parallel  to  the  c-axls  than  normal  to  it.  This  is  because 
brine  and  air  bubble  layers  are  regions  of  lower  conductivity,  and  more 
of  these  layers  are  Intersected  by  heat  conducted  along  the  c-axls  than 
perpendicular  to  it.  Since  horizontal  temperature  gradients  in  natural 
sea  ice  are  usually  negligible,  only  heat  flow  in  the  vertical  (i.e., 
perpendicular  to  the  c-axis)  need  be  considered  for  practical  purposes. 
Detailed  theories  for  k^,  taking  into  account  such  factors  as  brine 
pocket  distribution,  air  bubbles,  and  the  precipitation  of  solid  salts, 
have  been  worked  out  (Anderson,  1958;  Schwerdtfeger,  1963).  Tables  of 
ki(Ti,Si,Pi)  are  given  in  Schwerdtfeger  (1963).  An  approximate 
expression  for  k^  was  proposed  by  Untersteiner  (1961), 

BSi 

=  K  +  ^  ^  (16) 

^1 

where  Tj^  is  in  °C  and  6  -  3.11’10'^  cal  cm’l  s*l  (0.13  W  m'^)  is  a 
constant.  Although  somewhat  crude,  (16)  is  particularly  convenient  in 
applications  requiring  a  simple  functional  relationship  between  k^,  T^, 
and  S^.  Predicted  values  of  kj^(S^,T£)  are  shown  in  Fig.  12.  Note  that 
there  is  little  dependence  on  temperature  below  about  -4°C.  Under 
natural  conditions,  k^  may  vary  by  10-20%  in  perennial  ice  and  by  up  to 
50%  in  young  ice. 

At  first  glance  it  would  appear  that  (16)  predicts  an  infinite 
value  for  k^  as  Tj^  -►  0°C.  This  turns  out  not  to  be  the  case  because  Tj^ 
cannot  approach  0°C  unless  ->  0.  Whereas  pure  ice  has  a  well  defined 
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Fig.  12.  Dependence  of  thermal  conductivity  in  sea  ice  on  ice 
salinity  (°/qq)  and  temperature  according  to 
Untersteiner  (1961). 

temperature  (0°C)  at  which  phase  transitions  occur,  sea  ice  does  not. 
Any  change  in  produces  internal  phase  changes.  It  can  be  shown, 
however,  that  for  ice  of  a  given  salinity  there  is  a  temperature  (T|g) 
above  which  the  ice  cannot  exist  without  completely  melting  away. 
According  to  Ono  (1967),  =  -0.05411  S^.  This  is  consistent  with  the 

notion  that  the  temperature  of  sea  ice  can  reach  0°C  only  if  -  0. 


2.  Specific  Heat 

The  specific  heat  of  a  substance  is  defined  as  the  amount  of  heat 
necessary  to  raise  a  unit  mass  of  the  material  1°C,  and  Is  a  fundamental 
property  needed  In  calculating  time -dependent  temperature  changes.  The 
specific  heat  of  sea  Ice  (c^)  Is  the  sum  of  the  heat  required  to  warm 
the  Ice  and  the  brine  plus  that  needed  for  melting  In  the  brine  pockets. 
Ono  (1967)  has  derived  a  theoretical  expression  for  c^  applicable  to 
temperatures  above  -8°C: 

bSi 

=  Co  +  aT^  +  =7-2-  ,  (17) 

where  Cq  “  0.505  cal  g*^  (2113  J  kg*^  ®C*^)  is  the  specific  heat  of 

pure  Ice,  a  -  0.0018  cal  g*^  ®C*^  (7.53  J  kg‘^  ®C*^)  is  a  constant,  and 
b  -  4.3115  cal  ®C  g*^  (0.018  MJ  kg’^)  is  also  a  constant.  A  slightly 
simpler  form  which  gives  similar  results  is  the  empirical  formula  of 
Vntersteiner  (1961):  “  c©  where  ^  -  4.1  cal  ®C  g'^ 

(0.0172  MJ  ®C  kg*^).  Tables  and  analytical  expressions  are  also 
available  (Anderson,  1958;  Scbwerdtfeger,  1963)  for  colder  Ice.  Values 
of  c^(T£,S^)  are  shown  in  Fig.  13.  As  with  k^,  the  largest  changes  In 
Cj^  occur  at  temperatures  above  -4  to  -6°C,  but  the  magnitude  of  the 
change  Is  much  larger.  Under  natural  conditions  c^  can  vary  by  a  factor 
of  20  or  more.  The  figure  demonstrates  that  It  becomes  Increasingly 
difficult  to  change  T^  as  T^  ->■  primarily  because  mass  changes  In  the 
brine  pockets  consume  an  Increasing  share  of  the  available  energy. 

3.  Latent  Heat  of  Fusion 

The  latent  heat  of  fusion  (L^)  of  pure  ice  is  the  amount  of  heat 
required  to  melt  (or  freeze)  a  unit  mass  of  Ice  at  0°C.  Because 
Internal  melting  and  freezing  occur  continuously  over  a  wide  temperature 
range,  the  traditional  concept  of  absorption  or  release  of  latent  heat 
at  a  constant  temperature  does  not  strictly  apply  to  sea  Ice. 
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where  -  79.68  cal  g  (0.3335  MJ  kg'^).  Simpler  expressions  are  also 
possible.  Since  the  amount  of  ice  present  in  a  unit  mass  of  sea  ice  is 
(1  -  O.Olv^),  -  (1  -  0.01v^)Lq,  and  from  Eq.  (13) 


[1  -  8^(0.000532  - 


These  equations  can  be  used  to  estimate  the  heat  of  formation  of  sea 
ice.  Suppose,  for  example,  that  salinity  near  the  bottom  of  a  growing 
ice  cover  is  found  to  be  8°/oo.  Then,  if  Tf  -  -1.8°C,  s  22%  and  - 

61.2  cal  g  where  represents  the  amount  of  heat  that  must  be 
conducted  away  from  the  growing  interface  to  form  one  gram  of  salty  ice . 
As  the  ice  cools  below  T£,  additional  latent  heat  is  released  by 
freezing  in  the  brine  pockets,  but  that  heat  is  accounted  for  in  c^  and 
is  not  included  in  the  heat  of  formation. 


4.  Diffusivity  and  Conductance 

A  parameter  appearing  frequently  in  the  literature  is  the  thermal 
diffusivity ,  -  ^i/*^i®i’  ®  simple  combination  of  previously  discussed 

properties.  This  combination  arises  naturally  from  the  form  of  the  heat 
conduction  equation  (see  Section  5B)  and  can  be  calculated  directly  from 
observations  of  the  rate  of  change  of  temperature  profiles  in  the  ice 
(Schwerdtfeger ,  1963).  There  is  a  very  strong  decrease  in  with 
temperature  above  -4®C.  The  most  complete  study  of  k:£(T£,S£)  for  sea 
ice  was  made  by  Ono  (1968). 

A  useful  quantity  for  estimating  growth  rates  and  heat  conduction 
in  snow-covered  first-year  ice  is  the  thermal  conductance,  y  -  h^kg/ 
(k^hg  +  kgh) .  Because  kg  «  k^,  small  amounts  of  snow  can  have  a  large 
effect  on  the  amount  of  heat  conducted  to  the  surface  (F^.).  In 
situations  where  hg  undergoes  significant  changes  in  time  or  space,  the 
insulating  effect  of  the  snow  cannot  be  ignored.  Equation  (5)  shows 
that  surface  temperature  and  y  are  sufficient  to  obtain  F^,  when  the  ice 
is  cold.  If  F^  is  known  or  assumed  to  be  negligible,  the  amount  of 


growth  can  be  calculated  directly  from  F^.  and  L^.  The  dependence  of  Y 
on  h  and  hg  is  shown  in  Fig.  14  for  values  of  and  kg  typical  of 
winter  conditions  in  the  Arctic. 


B.  Mechanical  Properties 

Sea  ice  is  subject  to  a  variety  of  forces  arising  from  the  movement 
of  air  and  water  across  its  horizontal  boundaries ,  and  from  Internal 
stresses  generated  by  the  surrounding  ice.  The  resulting  deformations 
are  of  two  types:  elastic  strains  which  are  reversible  when  the  stress 
is  removed  and  inelastic  strains  which  are  irreversible.  Ice  is  often 
described  as  a  visco-elastic  material,  that  is,  when  deformed  rapidly  or 
for  brief  periods  it  behaves  elastically,  while  for  sustained  loadings 
it  tends  to  creep  with  a  rate  depending  on  time  and  the  amount  of 
stress.  The  mechanical  behavior  on  nonsaline  ice  depends  on  tempera¬ 
ture,  porosity,  grain  size,  strain  rate,  and  structure.  Sea  ice  is 
further  complicated  by  the  presence  of  brine  Inclusions,  strong  crystal 
orientation,  biological  activity,  and  organized  pore  geometry. 

Much  of  what  is  known  about  the  rheology  and  strength  of  ice  has 
been  learned  from  laboratory  measurements.  Continual  variations  in  the 
state  of  the  material  make  such  measurements  difficult,  and  much  work 
remains  to  be  done.  It  is  not  the  intent  here  to  summarize  the  numerous 
(and  sometimes  contradictory)  experimental  data.  Rather,  we  shall  try 
to  present  a  more  qualitative  picture  of  how  sea  ice  responds  to 
various  types  of  mechanical  forcing,  and  how  this  response  depends  on 
other  physical  properties  of  the  ice.  An  excellent  review  of  current 
knowledge  regarding  the  mechanical  behavior  of  sea  ice,  together  with  a 
discussion  of  the  underlying  continuum  mechanics  concepts,  has  been 
given  by  Mellor  (1983).  Additional  information  on  the  small  scale 
physics  can  be  obtained  from  Doronin  and  JCheisin  (1975). 

1 .  Elastic  Constants 

The  theory  of  elasticity  forms  the  basis  for  studies  of  deformable 
solids.  A  perfectly  (or  linearly)  elastic  material  is  one  in  which 
strain  (c)  in  the  direction  of  the  applied  stress  (a)  is  directly 


Fig.  14. 


ICE  THICKNESS  (m) 

Thermal  conductance  (W  m*^  ®C*^)  of  the  arctic  ice  cover 
during  winter  and  spring  as  a  function  of  ice  and  snow 
thickness  (after  /faykut,  1978). 


proportional  to  a,  i.e.,  a  -  Ee,  where  E  h  Young's  modulus,  the  ratio  of 
axial  stress  to  strain.  The  ratio  of  the  lateral  strain  to  the  axial 
strain  is  called  the  Poisson  ratio  (y) .  Together,  E  and  y  characterize 
the  elastic  behavior  of  a  material.  In  nonsaline  ice,  E  depends  only 
weakly  on  Tj^,  but  porosity  (air  bubble  volume)  has  a  significant  effect. 
In  sea  ice,  Tj^  strongly  affects  the  porosity  (brine  volume)  and  hence  E. 
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Small  scale  measurements  (Langleben  and  Pounder,  1963)  on  first-year  ice 
over  a  brine  voliame  range  of  0-8.5%  indicate  a  relationship  of  the  form, 
E  -  a  +  bv^j,  where  a  and  b  are  constants.  These  measurements  also 
suggest  that  p  is  relatively  constant  ('^<0.3)  in  sea  ice  and  Independent 
of  and  Theoretical  studies  of  the  dependence  of  E  on  have 

been  carried  out  by  Bergdahl  (1977)  assuming  that  brine  pockets  and 
channels  within  the  ice  could  be  modeled  as  vertical  rectangular  tubes . 
For  stresses  along  the  tubes  he  found  a  relationship  with  the  same  form 
as  the  empirical  studies.  However,  for  stresses  normal  to  the  tubes, 
the  relationship  was  of  the  form  E  -  a/(l  -  bv^^) ,  where  b  depends  on  the 
assumed  geometry  of  the  brine  pockets .  This  result  is  supported  by 
numerous  observations  indicating  that  the  direction  of  applied  stress  is 
an  Important  factor  in  the  mechanical  response  of  natural  sea  ice. 

Conditions  under  which  the  ice  behaves  elastically  are  fairly 
limited  and  often  are  not  satisfied  in  field  and  laboratory  measurements 
reported  in  the  literature.  Nevertheless,  elastic  analysis  is  still 
applied  in  many  problems  that  involve  inelastic  processes.  Elastic 
moduli  determined  under  these  conditions  are  called  "effective"  moduli 
and  are  denoted  by  E'  and  u'.  Effective  moduli  have  proved  to  be  useful 
for  many  practical  purposes,  and  a  variety  of  studies  have  been  made  to 
determine  their  dependence  on  T^,  porosity,  strain  rate,  and  stress 
state.  In  general  E'  is  substantially  less  than  E,  except  at  very  low 
temperatures  and  high  rates  of  stress  or  strain.  Increasing  porosity 
causes  a  much  more  rapid  decline  in  E'  than  in  E,  there  being  roughly  a 
50%  decrease  in  E'  as  temperature  increases  from  -30  to  -5°C.  Values 
of  p'  at  very  low  stress  rates  approach  0.5  and  vary  between  0.35-0.4 
for  high  stress  rates.  In  contrast  to  fresh  ice,  there  is  a  tendency  in 
saline  ice  for  p'  to  Increase  with  increasing  T^  (and  hence  .  In 
natural  sea  ice,  where  observations  {Wang,  1981)  indicate  that  p'  equals 
0-0.2  in  the  vertical  and  0.8 -1.2  in  the  horizontal,  sample  orientation 
is  of  utmost  importance. 
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2.  Strength 

The  most  useful  property  for  engineering  purposes  is  strength,  the 
ratio  of  the  minimum  load  necessary  to  fracture  the  ice  to  its  cross- 
sectional  area.  Depending  on  temperature,  structure,  and  the  rate  of 
deformation,  either  brittle  or  plastic  fracture  may  occur.  Several 
different  parameters  are  commonly  used  to  characterize  ice  strength. 
Uniaxial  compressive  strength  (a^) ,  the  maximum  stress  that  the  ice  can 
sustain  in  a  particular  direction  without  fracturing,  is  usually 
obtained  by  compression  of  ice  cylinders  at  a  specific  rate  of  strain. 
Results  show  that  in  sea  ice  increases  as  strain  rate  (e:)  to  approxi¬ 
mately  the  one -fourth  power.  As  expected,  decreases  more  strongly 
with  temperature  in  sea  ice  than  in  either  lake  or  river  ice ,  and  is 
much  larger  when  force  is  applied  in  the  direction  of  ice  growth 
(vertical)  than  when  applied  in  the  horizontal.  The  ice  is  weakest  when 
the  stress  is  45°  to  the  horizontal  plane.  Uniaxial  tensile  strength 
(o^)  can  be  measured  by  applying  tension  to  carefully  designed  ice 
specimens,  but  practical  difficulties  with  this  approach  have  led  most 
Investigators  to  infer  from  observations  of  beam  flexure.  The 
behavior  of  both  small  ice  beams  in  the  laboratory  and  large  ice  canti¬ 
levers  sawed  in  natural  sea  ice  sheets  has  been  studied.  At  low  strain 
rates  (  e  <  10  s  ) ,  a.j,  and  0^  are  roughly  equal;  however,  for  high 
rates  of  loading  Oj.  is  much  smaller  than  As  with  increases 

with  decreasing  temperature  and  salinity  and  is  considerably  greater  for 
vertically  oriented  specimens  of  columnar  ice  than  for  horizontal  speci¬ 
mens.  Other  measures  of  ice  strength  (e.g.,  flexural  strength,  shear 
strength,  fracture  toughness,  bearing  capacity)  have  also  been  investi¬ 
gated.  Results  from  much  of  this  work  have  been  siommarlzed  by  Doronin 
and  Kheisin  (1975)  and  Mellor  (1983).  Little  has  been  published  on  the 
mechanical  properties  of  multiyear  ice.  What  evidence  is  available 
indicates  that  its  strength  is  highly  variable,  with  differences  of  an 
order  of  magnitude  being  observed  in  ice  with  similar  strain  rates. 

Much  more  work  is  needed  before  it  will  be  possible  to  draw  comparisons 
with  the  results  from  first-year  ice. 


Brine  pockets  and  channels  weaken  sea  ice,  and  fracturing  tends  to 
follow  planes  that  pass  through  these  features.  This  tendency  has  been 
used  to  simplify  the  treatment  of  pore  geometry  and  fracture  strength  in 
theoretical  sea  ice  models  (Weeks  and  Assur,  1969;  Weeks  and  Ackley, 
1982).  In  these  models  brine  pockets  are  assumed  to  be  vertically 
oriented  elliptical  cylinders,  elongated  in  the  direction  of  the  plate¬ 
lets,  which  occur  in  parallel  rows.  Model  results  show  a  relationship 
between  fracture  strength  and  of  the  form  -  ^^^(1  -  aVj^^)  ,  where 
is  the  fracture  strength  of  sea  ice  containing  no  brine ,  and  a  and  k 
depend  on  the  arrangement  of  the  brine  pockets  and  on  how  their  shape 
changes  with  v^.  Assuming  that  pockets  enlarge  in  all  directions  to 
maintain  a  constant  geometry  results  in  k  -  2/3,  while  assuming  that 
expansion  of  the  pockets  primarily  affects  the  degree  of  elongation 
results  in  k  -  1.  If,  on  the  other  hand,  it  is  assumed  that  the  rela¬ 
tive  length  of  the  brine  columns  remains  constant  and  that  only  the 
cross-sectional  area  changes  with  to  maintain  the  same  geometrical 
shape,  it  follows  that  k  -  1/2,  in  agreement  with  observations  (Weeks 
and  Ackley,  1982). 

A  subject  that  has  received  surprisingly  little  attention  is  the 
effect  of  biological  activity  on  the  strength  and  breakup  of  an  ice 
cover.  It  has  long  been  known  that  many  species  of  microalgae  Inhabit 
sea  ice,  but  the  degree  to  which  their  presence  may  alter  the  mechanical 
properties  of  the  ice  does  not  appear  to  have  been  generally  appreci¬ 
ated.  A  particularly  Interesting  series  of  observations  was  carried  out 
by  Buynitskiy  (1968)  during  summer  cruises  in  the  Antarctic.  He  found 
that  substantial  quantities  of  algae  were  present  in  essentially  all  the 
ice  between  the  extreme  edge  and  the  shore.  Cell  densities  as  high  as 
33,000  per  cubic  centimeter  were  recorded,  with  the  average  being  about 
4300  per  cubic  centimeter,  roughly  an  order  of  magnitude  larger  than  in 
the  underlying  water  column.  Late  in  the  summer  diatoms  occurred 
throughout  the  ice  cover,  while  earlier  in  the  year  they  tended  to 
concentrate  in  brine  channels  and  in  horizontal  layers  ranging  from  a 
few  centimeters  to  tens  of  centimeters  in  thickness.  The  dark  color  of 
the  algae  enhances  absorption  of  shortwave  radiation,  accelerating  the 


rate  of  Internal  melting  near  the  algae  and  thereby  weakening  the  Ice. 
Ship  travel  through  such  Ice  was  found  to  be  substantially  easier  than 
through  Ice  with  little  algae.  In  contrast  to  pure  sea  Ice  which  splits 
mainly  In  directions  normal  to  the  plane  of  freezing,  antarctic  Ice  was 
observed  to  cleave  frequently  through  horizontal  planes  populated  by  the 
algae.  Strength  tests  were  carried  out  on  various  Ice  samples  contain¬ 
ing  algae.  It  was  found  that  the  flexural  strength  of  such  Ice  was 
reduced  by  an  average  of  40%  over  that  of  pure  sea  Ice ,  with  the 
greatest  reduction  In  strength  occurring  for  stresses  applied  perpen¬ 
dicular  to  the  freezing  Interface.  The  bearing  capacity  of  this  Ice  was 
thus  seriously  decreased.  Compressive  tests  Indicated  that  o^  was 
almost  unchanged  by  the  presence  of  diatoms  for  vertically  applied 
stresses,  but  decreased  by  a  factor  of  2  for  horizontal  stresses. 
Information  In  the  published  account  Is  Inadequate  to  Identify  the  exact 
mechanisms  whereby  the  algae  achieve  this  weakening,  but  It  Is  probably 
due  more  to  warmer  temperatures  and  higher  brine  volumes  In  the  vicinity 
of  the  horizontal  layers  than  to  any  direct  effect  on  the  microstructure 
of  the  Ice.  Whether  biological  activity  Is  high  enough  anywhere  In  the 
Arctic  to  produce  similar  effects  Is  uncertain. 

3 .  Large  Scale  Behavior 

Although  the  large-scale  mechanical  behavior  of  a  sea  Ice  cover  can 
be  vastly  different  from  that  observed  In  small  laboratory  samples,  the 
latter  provide  a  basis  for  understanding  and  treating  large-scale 
behavior.  Small-scale  strength  measurements,  for  example,  Indicate  that 
tensile  forces  generated  by  wind  and  water  stresses  are  much  too  small 
to  crack  a  floating  sheet  of  unflawed  Ice.  Fracture  analysis  similar  to 
that  used  with  the  laboratory  data  shows  that  leads  In  natural  sea  Ice 
must  form  as  a  result  of  the  growth  of  small  existing  cracks.  Such 
cracks  can  be  produced  by  thermal  stresses  or  wave  action  and  are  very 
common  In  natural  Ice  covers.  Compressive  stresses  are  likewise  Inade¬ 
quate  to  cause  wide-scale  crushing  In  an  Ice  sheet,  although  they  may 
produce  buckling  In  areas  of  very  thin  Ice.  This  suggests  that  other 
processes  must  be  Involved  In  the  formation  of  pressure  ridges.  A 
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likely  factor  contributing  to  ridge  building  is  overthrusting  of  the  ice 
along  a  discontinuity  (such  as  an  old  lead),  causing  flexural  breakage 
to  occur  at  much  lower  stress  levels  than  needed  for  buckling.  Over¬ 
thrusting  should  take  place  primarily  in  thinner  ice  and  does  not 
explain  how  thicker  ice  can  become  involved  in  ridging.  A  possible 
explanation  for  this  can  be  found  in  a  ridging  model  developed  by 
Parmerter  and  Coon  (1972) .  Thin  ice  in  an  existing  lead  is  assumed  to 
break  up  as  the  lead  closes  in  response  to  compresslonal  forces.  The 
resulting  rubble  is  piled  on  and  beneath  the  surrounding  thick  ice, 
disturbing  its  isostatic  balance  which  then  leads  to  flexural  breakage 
adjacent  to  the  lead.  Incorporation  of  fragments  from  the  parent  ice 
into  the  rubble  pile  and  further  closure  allow  the  process  to  continue 
until  a  limit  height  is  attained.  Additional  activity  then  acts  to 
widen  the  ridge,  but  not  to  Increase  its  height.  Whether  this  repre¬ 
sents  a  quantitatively  accurate  simulation  of  the  ridging  process  is  not 
known;  however,  predicted  limit  values  appear  to  be  well  correlated  with 
ridge  height  data. 

Central  to  the  development  of  large-scale  models  of  ice  drift  and 
deformation  are  (1)  a  constitutive  law  relating  ice  stress  to  deforma¬ 
tion  and  (il)  some  measure  of  the  aggregate  strength  of  the  ice. 

Numerous  attempts  have  been  made  to  model  the  ice  as  a  linear  viscous 
fluid,  but  results  from  these  models  Indicate  that  nonlinear  effects  are 
important.  Physical  considerations  suggest  that  a  plastic  rheology 
should  provide  a  more  realistic  description  of  large-scale  dynamic 
behavior,  and  both  elastic-plastic  and  viscous -plastic  approaches  have 
been  used  with  some  success  (Coon  et  al,,  1974;  Hibler,  1979). 

The  overall  compressive  strength  (o^.  )  of  the  ice  pack  is  strongly 
affected  by  local  variations  in  thickness.  Converging  motions  cause  the 
thinnest  (weakest)  ice  to  ridge.  Increasing  its  resistance  to  further 
compression.  As  convergence  continues,  Increases  and  thicker  ice 
becomes  Involved  in  the  deformation.  The  magnitude  of  therefore 
depends  on  how  thinner  ice  is  distributed  within  the  region.  A  formal 
theory  has  been  developed  {Rothrock,  1975)  for  the  plastic  rheology  that 
relates  the  rate  of  work  done  on  the  ice  through  ridging  to  the  strain 


and  state  of  stress  providing  a  method  of  calculating  a^*  for  a  given 
distribution  of  ice  thickness.  Numerical  experiments  indicate  that 
values  of  are  typically  on  the  order  of  10^-10^  N  m”'*’.  A  more 
detailed  summary  of  work  in  this  area  can  be  found  in  Hibler  (1980) . 

C.  Optical  Properties 

The  optical  properties  of  a  material  define  how  it  interacts  with 
solar  radiation.  Shortwave  energy  exchange  over  the  polar  oceans  is 
complicated  by  nonuniformities  in  the  ice  cover  which  cause  its  optical 
properties  to  exhibit  large  variations  in  both  the  horizontal  and 
vertical.  At  the  height  of  the  melt  season,  for  example,  a  relatively 
small  area  may  contain  leads,  thin  ice,  surface  melt  ponds,  drained 
white  ice,  and  snow  patches,  all  having  significantly  different  optical 
properties.  Both  the  surface  and  internal  structure  of  the  ice  are 
altered  by  the  absorption  of  shortwave  radiation,  so  that  its  optical 
properties  undergo  temporal  variations  as  well  as  spatial.  Although  the 
optical  properties  of  the  ice  pack  do  change  gradually  during  the  winter 
months,  the  largest  changes  occur  during  the  summer  melt  cycle  with  the 
decay  of  the  snow  cover  and  the  formation  of  melt  ponds.  In  contrast  to 
snow  and  bare  ice,  ponds  are  areas  of  low  reflectivity  where  large 
amounts  of  energy  are  absorbed  and  transmitted  to  the  underlying  ice  and 
ocean.  Accompanying  the  surface  changes  are  changes  in  the  transparency 
of  the  ice  caused  by  Internal  melting  and  Increasing  brine  voluune. 
Quantitative  knowledge  of  how  light  absorption,  transmission,  and 
reflection  vary  in  time  and  space  is  of  fundamental  Importance  in  a  wide 
variety  of  studies  involving  the  ice  and  upper  ocean. 


1 .  Albedo 

The  albedo  (a)  is  defined  as  the  ratio  of  reflected  to  incident 
shortwave  radiation  and  is  a  measure  of  net  solar  input  at  the  surface. 
Total  (Integrated  over  all  wavelengths)  albedo  measurements  of  polar  ice 
and  snow  have  been  made  by  many  investigators  (see  Grenfell  and  Maykut , 
1977).  Albedos  of  cold  polar  snow  are  generally  in  the  0.80-0.85  range, 
dropping  to  0.70-0.75  at  the  onset  of  melting.  Snow-free  multiyear  ice 


has  an  albedo  of  0.70-0.75  when  the  surface  temperature  (T^)  is  below 
freezing  and  0.55-0.70  when  -  0°C.  Thicker  (1-2  m)  first-year  ice 
has  a  values  of  0.5-0. 6  when  cold  and  0.3-0. 5  when  melting.  The  albedo 
of  snow- free  young  ice  exhibits  a  strong  dependence  on  thickness  which 
may  be  approximated  by  a  -  0.44  h  +  0,08,  where  h  is  in  meters 
(Maykut ,  1982).  The  higher  reflectivity  of  multiyear  ice  is  due  to  the 
presence  of  a  decomposed  surface  scattering  layer  (typically  5-15  cm 
thick)  and  to  greater  vapor  bubble  density  in  the  underlying  ice .  Melt 
pond  albedos  depend  on  water  depth  and  the  state  of  the  underlying  ice . 
Deep  ponds  remaining  at  the  end  of  the  melt  season  can  have  albedos  as 
low  as  C.15,  while  values  of  0.4  are  typical  of  the  shallow  ponds  that 
form  early  in  the  melt  season.  Contamination  by  dirt  or  algae  can 
reduce  these  values  somewhat.  Lead  albedos  vary  with  wave  activity  and 
sun  angle ,  but  average  about  0.08. 

While  total  albedos  are  easy  to  measure  and  apply,  their  usefulness 
is  limited.  For  example,  analysis  of  satellite  observations,  estimation 
of  cloud  effects  on  a,  and  calculations  of  light  absorption  in  the  ice 
require  an  understanding  of  how  the  radiative  properties  of  the  ice  vary 
with  wavelength  (X).  Spectral  albedos  (a^)  have  been  measured  over  many 
surface  types  in  the  Arctic  (Grenfell  and  Maykut,  1977).  Values  for 
first-year  and  multiyear  ice  are  shown  in  Fig.  15.  Albedos  are 
relatively  constant  in  the  visible  region  (400-700  nm) ,  but  begin  to 
decrease  in  the  near  infrared  at  a  rate  that  is  related  to  the  liquid 
water  content  in  the  surface  layers.  Melting  first-year  ice  can  be 
divided  into  two  main  types:  (1)  water- saturated  blue  ice  (which 
Includes  shallow  melt  ponds),  and  (ii)  white  ice  consisting  of  a  5-10  cm 
thick  layer  of  drained  ice  underlain  by  blue  ice.  At  longer  wavelengths 
light  penetration  is  small  so  that  is  determined  by  the  properties  of 
the  surface  layer;  at  shorter  wavelengths  backscattering  from  ice 
beneath  the  surface  becomes  important.  This  explains  the  similarity 
between  in  first-year  white  ice  and  melting  multiyear  ice  at  longer 
wavelengths,  and  the  difference  at  shorter  wavelengths  that  results  from 
lower  backscattering  in  the  underlying  first-year  ice.  The  albedo  of 
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Fig.  15.  Spectral  albedos  observed  over  bare  sea  Ice  In  the 

Arctic:  (a)  frozen  multiyear  Ice,  (b)  melting  ualtlyear 

Ice,  (c)  melting  first-year  white  Ice,  and  (d)  melting 
first-year  blue  Ice  (after  Grenfell  and  Maykut,  1977). 

melting  blue  Ice  Is  substantially  lower  than  that  of  first-year  white 

Ice  at  all  wavelengths .  At  shorter  wavelengths  this  Is  due  to  the 

absence  of  the  white  scattering  layer,  and  at  longer  wavelengths  to 

Increased  liquid  water  In  the  surface  layers .  Spectral  albedos  for  snow 

and  melt  ponds  are  plotted  In  Fig.  16.  The  albedo  of  cold  dry  snow  Is 

high  and  has  only  a  weak  dependence  on  A.  Increasing  density  and  liquid 
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Fig.  16.  Spectral  albedos  observed  over  snow  and  melt  ponds  in 

the  Arctic:  (a)  dry  snow,  (b)  wet  new  snow,  (c)  melting 
old  snow,  (d)  frozen  melt  pond,  (e)  early  season  melt 
pond,  (£)  mature  melt  pond  on  multiyear  ice,  (g)  melt 
pond  on  £lrst>year  ice,  and  (h)  old  melt  pond  on 
multiyear  ice  (after  Grenfell  and  Haykut,  1977). 


water  content  in  the  snow  decrease  a^,  particularly  in  the  infrared. 
Pond  albedos  tend  to  be  high  at  shorter  wavelengths  with  a  sharp 
decrease  between  500  and  800  nm.  Because  water  in  the  ponds  is  rela¬ 
tively  transparent  at  short  wavelengths,  albedos  below  500  nm  are 
determined  primarily  by  the  scattering  properties  of  the  underlying  ice 
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Between  500  and  800  nm  becomes  increasingly  insensitive  to  the  under¬ 
lying  ice  as  absorption  in  the  water  becomes  the  dominant  factor.  Above 
800  nm  is  determined  only  by  Fresnel  reflection  from  the  pond 
surface.  Note  that  most  of  the  albedo  differences  between  individual 
ponds  occur  in  the  visible  region  and  are  readily  apparent  to  the  naked 
eye.  Ice  and  snow  albedos  decrease  strongly  above  1000  nm,  reaching 
values  of  0.1  or  less  at  1400  nm  after  which  the  decrease  becomes  much 
more  gradual  (Grenfell  and  Perovich,  1984). 

The  total  albedo  can  be  expressed  in  terms  of  wavelength-dependent 
quantities  as  a  =  (A)d^/(j  Fj. (A)dX,  where  Fj.(X)  is  the  downward 

flux  of  solar  radiation  reaching  the  surface  at  wavelength  X.  The 
presence  of  clouds  causes  measured  values  of  a  to  Increase.  The  reason 
is  that  clouds  remove  proportionately  more  energy  from  Fy(X)  at  long 
wavelengths  than  at  short  ones,  so  that  contributions  from  the  more 
reflective  shorter  wavelength  region  are  weighted  more  heavily  in  the 
Integral  when  the  sky  is  overcast.  Typical  increases  observed  for  the 
relatively  thin  clouds  found  in  the  Arctic  are  about  15%  (Grenfell  and 
Perovich,  1984).  Larger  increases  can  be  expected  in  regions  where 
clouds  have  greater  optical  depths. 

2.  Extinction  Coefficient 

Because  ice  is  a  translucent  material,  light  penetrates  into  the 
interior  where,  at  any  given  level,  some  is  absorbed  and  some  trans¬ 
mitted.  Light  absorbed  within  the  ice  produces  changes  in  the  brine 
volume  and  provides  energy  for  biological  activity.  Shortwave  radiation 
transmitted  completely  through  the  ice  is  largely  absorbed  in  the  upper 
10-20  m  of  the  ocean,  making  an  important  contribution  to  the  oceanic 
heat  flux  at  the  bottom  of  the  ice.  Light  attenuation  in  the  ice  is 
characterized  by  its  spectral  extinction  coefficient  (k^^)  defined  as 


F(X,z)  dz  ’ 


relatively  constant  values  of  in  the  400-500  nm  region,  followed  by 
strongly  increased  attenuation  at  longer  wavelengths.  The  greatest 
attenuation  was  observed  in  cold  dry  snow  where  were  about  20  times 

as  large  as  those  within  the  ice.  The  "transition  zone"  in  Fig.  17  is 
a  region  of  rapidly  changing  Kx  that  occurs  within  the  5-15  cm  thick 
surface  scattering  layer  on  multiyear  ice  as  it  grades  into  the  more 
uniform  interior  ice.  Differences  between  '^x'^  first-year  and 
multiyear  ice  result  from  differences  in  brine  volume  and  bubble 
density.  Minlmvim  values  (curve  f)  were  found  beneath  a  mature  melt  pond 
where  brine  volume  in  the  upper  part  of  the  ice  exceeded  50%. 

The  strong  attenuation  at  longer  wavelengths  causes  the  spectral 
composition  of  the  light  to  undergo  rapid  changes  as  it  passes  through 
the  ice.  Measurements  of  downwelllng  Irradiance  taken  beneath  first- 
year  sea  ice  (Maykut  and  Grenfell ,  1975)  show  energy  maxima  occurring  in 
the  470-480  nm  wavelength  band,  corresponding  to  the  peak  in  the  inci¬ 
dent  solar  spectrum.  Red  wavelengths  are  quickly  absorbed  in  the 
surface  layers  of  the  ice  so  that  the  amount  of  energy  reaching  the 
ocean  is  negligible  above  700-800  nm,  even  when  the  ice  is  very  thin. 
Maximum  transmission  by  the  ice  occurs  in  the  450-550  nm  region,  regard¬ 
less  of  the  surface  conditions  or  ice  thickness .  The  nature  of  the 
surface  is,  however,  the  dominant  factor  in  determining  the  amount  of 
energy  reaching  the  ocean.  These  points  are  illustrated  by  Figs.  18 
and  19  which  show  the  spectral  distribution  of  light  beneath  different 
thicknesses  of  first-year  melt  pond/blue  ice  and  white  ice.  Because  the 
spectral  composition  of  the  Incident  radiation  field  is  rarely  measured, 
it  is  useful  to  be  able  to  estimate  F(X,h)  from  more  commonly  available 
data.  Values  in  Figs.  18  and  19  have  thus  been  expressed  as  a  percent 
of  the  total  Incoming  shortwave  radiation  measured  under  typical  arctic 
cloud  conditions .  The  presence  of  snow  on  the  ice  drastically  reduces 
the  amount  of  transmitted  light  as  Illustrated  in  Fig.  20. 

In  many  regions  a  thick  layer  of  algae  forms  on  the  underside  of 
the  ice  in  late  spring,  intercepting  most  of  the  transmitted  light.  With 
the  onset  of  bottom  ablation,  this  layer  tends  to  slough  off;  however. 


Wovelength  (nm) 

Fig.  18.  Spectral  distribution  of  downwelling  irradiance  beneath 
first-year  melt  pond/blue  ice  (after  MayJcut  and 
Grenfell,  1975).  Isopleths  are  expressed  in  percent  per 
nanometer  of  the  total  Incident  irradiance  at  the 
surface . 

algae  remaining  within  the  ice  are  still  sufficient  to  affect  the  spec¬ 
tral  distribution  of  the  transmitted  light.  The  characteristic  signa¬ 
ture  of  these  algae  is  a  broad  peak  in  the  transmission  spectra  between 
450-550  nm  and  a  secondary  peak  or  shoulder  at  about  540  nm  (Haykut  and 
Grenfell ,  1975).  Bottom  melting  and  brine  drainage  act  to  reduce  the 
density  of  interstitial  algae,  and  the  transmission  spectra  slowly 
return  to  a  fairly  narrow  maximum  at  470-480  nm  later  in  the  melt 
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Fig.  19.  Spectral  distribution  of  downwelllng  Irradlance  beneath 
first-year  white  ice  (after  Maykut  and  Grenfell,  1975). 

season.  Differences  in  transmission  spectra  measured  with  and  without 
algae  Indicate  that  the  algae  cause  maximum  absorption  at  435  nm,  with  a 
secondary  peak  at  about  670  nm. 

The  spectral  distribution  of  shortwave  radiation  within  the  ice  can 
be  calculated  by  integrating  Eq.  (20)  to  obtain 
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Fig.  20.  Total  dovnwelllng  irradiance  beneath  first-year  arctic 
sea  ice  as  a  function  of  ice  thickness  and  snow  depth 
(after  Maykut  and  Grenfell,  1975).  Isopleth  values  are 
given  as  a  percentage  of  the  total  incident  irradiance 
at  the  surface  under  overcast  conditions. 


where  f(A,0)  -  (1  -  a;^)Fp(A)  -  the  net  surface  flux  of  shortwave  radi¬ 
ation  at  wavelength  A.  Equation  (21)  is  known  as  Beer's  law.  In  most 
practical  applications  it  is  the  total  flux,  =  /*F(A,z)dA,  which  is 
of  interest.  By  defining  a  bulk  extinction  coefficient  (k^)  as 
<2  *  “  /q *^^^dA//^  F(A,0)e  "^^^dA,  it  is  possible  to  rewrite 
Beer's  law  as 
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where  -  /QF(^,o)dX  -  (l-a^)Fj.(X)dX  is  the  net  total  shortwave 
radiation  at  the  surface.  Grenfell  and  Haykut  (1977)  have  examined  the 
variation  of  with  depth  in  both  first-year  and  multiyear  ice.  They 
found  that  decreases  by  as  much  as  a  factor  of  50  in  the  upper  10-20 
cm  of  the  ice,  after  which  it  exhibits  only  a  weak  dependence  on  z. 

Very  little  difference  was  found  between  in  the  Interior  of  first- 

year  and  multiyear  ice.  For  purposes  of  estimating  F2,  the  ice  can  be 
thought  of  as  consisting  of  two  layers:  an  upper  layer  that  filters 
out  nearly  all  the  energy  at  longer  wavelengths ,  and  a  lower  layer  where 
<2  does  not  depend  strongly  on  either  ice  type  or  cloud  conditions. 
Shortwave  energy  absorbed  in  the  upper  layer  is  assumed  to  contribute 
immediately  to  temperature  or  mass  changes  at  the  surface,  while  energy 
absorbed  in  the  lower  layer  returns  slowly  to  the  surface  through 
conduction.  This  picture  provides  a  simple  method  of  estimating  F2.  If 
we  define  to  be  the  fraction  of  F^  that  passes  through  the  upper 
layer ,  then 
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where  Ic  is  the  average  bulk  extinction  coefficient  in  the  lower  layer 
and  Zq  is  the  thickness  of  the  upper  layer.  Equation  (23)  provides  the 
best  fit  to  observations  when  F  -  0.015  cm'^.  The  value  of  depends 
on  both  the  ice  type  and  the  spectral  distribution  of  the  incident 
radiation  (i.e.,  on  cloudiness).  Using  z^  -  10  cm,  Grenfell  and  Maykut 
(1977)  calculate  that  for  melting  white  ice  is  0.35  when  overcast  and 
0.18  when  clear;  for  blue  ice  it  is  0.63  when  overcast  and  0.43  when 
clear.  With  (23)  it  is  possible  to  calculate  the  amount  of  energy 
transmitted  to  the  ocean  by  setting  z  -  h,  while  the  amount  of  energy 
absorbed  at  a  depth  z(z  >  z^)  is  simply  kF^. 
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D.  Ice  Thickness  and  Floe  Size  Distributions 

It  was  pointed  out  earlier  that  continual  deformation  and  growth 
produce  an  ice  cover  that  Is  a  mixture  of  many  different  thicknesses  of 
Ice.  Because  many  of  the  physical  properties  of  the  ice  cover  depend 
strongly  on  h,  it  is  useful  to  have  a  way  to  quantify  thickness 
variations  within  a  particular  region.  To  accomplish  this,  Thorndike  et 
al.  (1975)  have  introduced  a  probability  density  function  g(h)  called 
the  ice  thickness  distribution,  where  g(h)dh  is  the  fractional  area 
covered  by  ice  in  the  thickness  range  hj^  <  h  <  h2.  Processes  that  act 
to  change  g(h)  are  thermodynamic  growth,  ice  advectlon,  ridging,  and 
divergence.  A  thickness  distribution  model  that  includes  these 
processes  has  been  developed  and  tested  using  data  from  the  Central 
Arctic  (.Thorndike  et  al.,  1975).  While  producing  reasonable  results, 
the  model  parameterizes  a  number  of  complex  and  poorly  understood 
processes  related  to  the  mechanical  redistribution  of  h.  Considerable 
work  on  the  mechanical  response  of  the  ice  to  various  strain  states 
needs  to  be  carried  out  before  it  will  be  possible  to  specify  these 
processes  with  any  degree  of  confidence.  Nevertheless,  g(h)  provides  a 
framework  for  coupling  the  dynamics  and  thermodynamics  of  the  ice  pack 
that  has  proved  to  be  very  useful  in  a  variety  of  modeling  efforts  and 
large-scale  studies  of  the  ice. 

Although  the  large-scale  behavior  of  the  ice  is  directly  influenced 
by  the  distribution  of  ice  thickness,  data  on  g(h)  are  difficult  to 
obtain.  At  present  the  only  practical  way  to  measure  g(h)  directly  is 
with  upward- looking  submarine  sonar.  Figure  21  shows  distributions  of 
ice  draft  derived  from  submarine  transects  north  of  Fram  Strait  in 
October  (Wadhams,  1981).  As  might  be  expected,  there  is  a  strong  peak 
at  about  3  cm,  close  to  the  thermodynamic  equilibrium  thickness,  and  a 
secondary  peak  at  about  50  cm  originating  from  young  ice  formed  earlier 
in  the  fall .  Level  (presumably  undeformed)  ice  makes  up  about  50%  of 
the  ice  in  the  region.  The  distribution  of  deformed  ice  follows  a 
negative  exponential  at  drafts  greater  than  8-10  m.  Measurements  in  the 
Beaufort  Sea  (Wadhsuns  and  Home,  1980)  show  similar  characteristics, 
although  the  secondary  peak  in  g(h)  occurs  at  about  2  m  rather  than  50 
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ICE  DRAFT  (m) 


October  ice  thickness  distributions  observed  in  the 
Arctic  Ocean  north  of  Fram  Strait  (after  Wadhams,  1981) 
(a)  the  distribution  of  level  ice,  and  (b)  the  complete 
distribution  Including  ridged  ice. 


cm,  probably  representing  second-year  ice  from  the  previous  season. 
While  a  systematic  analysis  of  spatial  variations  In  g(h)  has  not  be 
carried  out,  LeShack  (see  Hibler,  1980)  has  utilized  submarine  data 
calculate  the  average  thickness  (h)  at  many  locations  In  the  Arctic 
Basin  (Fig.  22).  The  dominant  feature  of  the  h  distribution  is  a  sh 
Increase  In  thickness  north  of  Greenland  and  the  Canadian  Archlpelag 
that  is  almost  certainly  the  result  of  ereater  rldeine  activity  In  t 
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While  providing  adequate  definition  of  variations  in  the  amount  of 
thicker  ice,  submarine  sonar  data  are  not  sufficiently  dense  in  time  to 
resolve  the  large  and  rapid  changes  that  occur  in  the  distribution  of 
young  ice.  Yet  knowledge  of  temporal  changes  in  g(h  <  1  m)  is  of 
fundamental  Importance  because  such  areas  not  only  make  a  major 
contribution  to  regional  ice  production  and  heat  exchange,  but  also 
control  the  mechanical  behavior  of  the  ice  pack.  Without  routine 
submarine  transects  or  satellite  Instrvunents  capable  of  resolving  young 
ice  thicknesses  to  within  10-20  cm,  g(h)  must  be  inferred  from  other 
data.  The  most  feasible  approach  is  to  predict  g(h)  theoretically  by 
combining  models  of  ice  growth  and  thickness  distribution  with  heat 
balance  and  ice  velocity  data.  Thermodynamic  forcing  for  the  models  is 
often  derived  from  climatology,  modified  on  the  basis  of  buoy  or 
satellite  observations.  Information  on  the  velocity  field  of  the  ice 
can  be  obtained  directly  from  the  motion  of  drifting  buoys  (Colony  and 
Thorndike,  1983),  or  calculated  from  atmospheric  pressure  data  using  a 
dynamic  ice  model  (Coon  et  al.,  1974;  Hibler,  1979). 

A  geometric  property  of  the  ice  cover  that  has  received  relatively 
little  attention  is  the  floe  size  distribution,  N(d) ,  defined  as  the 
number  of  floes  of  mean  diameter  d  per  unit  area.  Floe  size 
distribution  is  believed  to  be  a  significant  factor  in  problems  related 
to  the  summer  ice  pack.  Fart  of  the  shortwave  energy  deposited  in  the 
upper  few  meters  of  a  lead  goes  to  lateral  melting  on  floe  edges  causing 
a  decrease  in  ice  concentration  (A^)  and  an  Increase  in  the  amount  of 
open  water  (A^,) .  This  positive  feedback  between  solar  input  to  the 
ocean  and  the  rate  of  change  of  A^  can  have  a  large  effect  on  the  summer 
decay  and  retreat  of  the  ice  (Langleben,  1972).  As  the  ice  cover 
becomes  more  broken  up,  there  is  a  decreasing  probability  that  heat 
absorbed  in  the  leads  will  be  lost  to  the  atmosphere  or  underlying 
ocean,  and  a  corresponding  Increase  in  the  amount  of  heat  transferred  to 
floe  edges.  Simple  model  experiments  (Perovich,  1983)  show  that,  for  a 
given  ice  concentration,  the  rate  of  lateral  decay  depends  directly  on 
the  number  of  leads  present,  i.e.,  on  the  total  floe  perimeter.  Ice 
concentration  alone  is  thus  inadequate  to  characterize  the  state  of  the 


summer  ice  cover.  N(d)  provides  a  measure  of  total  floe  perimeter  and, 
unlike  g(h) ,  is  relatively  simple  to  obtain  from  aerial  photography. 
N(d)  appears  to  be  particularly  important  in  modeling  the  Marginal  Ice 
Zone,  where  it  affects  not  only  lateral  melting  but  also  drag  coeffi¬ 
cients,  melt  pond  coverage,  and  open  water  distribution.  The  equation 
governing  N(d)  is  analogous  to  that  for  g(h) ,  containing  terms  to 


describe  advection,  thermodynamics  (lateral  melting),  and  mechanical 
processes  (breaking) .  Except  near  the  ice  edge  where  wave  activity 
fractures  the  ice  in  a  predictable  fashion,  processes  involved  in  floe 
breakup  are  not  well  understood.  Measurements  of  floe  size  distribution 
have  been  made  in  several  locations  (Wadhams,  1973;  Vinje,  1977; 

Rothrock  and  Thorndike ,  1984).  Figure  23,  for  example,  shows  mid-August 
values  in  the  central  Beaufort  Sea.  For  this  case  the  decrease  in  N 
with  d  roughly  follows  a  power  law  between  0.05  and  10  km.  In  the 
Marginal  Ice  Zone,  N(d)  tends  to  assume  a  negative  exponential  distri¬ 
bution  (Vinje,  1977;  Wadhams,  1980),  but  in  the  interior  ice  N(d)  is 
much  more  variable  and  no  single  analytical  form  can  be  used  to  fit  all 
the  data  (Rothrock  and  Thorndike,  1984).  No  systematic  studies  have  yet 
been  carried  out  to  define  how  N(d)  evolves  during  the  summer  or  how 
this  evolution  differs  from  one  region  to  another. 

E.  Electromagnetic  Properties 

In  recent  years  there  has  been  Increasing  interest  in  the  electri¬ 
cal  properties  of  sea  ice  in  the  0.1-100  GHz  region.  Because  of  the 
extensive  cloud  cover  and  long  periods  of  darkness,  satellite  imagery  in 
the  visible  and  Infrared  are  of  limited  value  in  monitoring  the  state  of 
the  polar  ice  pack.  Signals  in  the  microwave  region,  however,  are  not 
affected  by  darkness  and  are  attenuated  only  slightly  by  clouds. 
Microwave  imagery  thus  provides  a  reliable  method  of  observing  the  ice 
throughout  the  year.  Both  active  and  passive  microwave  measurements 
have  been  carried  out  over  the  ice.  Active  methods  transmit  a  radar 
signal  to  the  ice  and  measure  the  resulting  reflection.  At  lower 
frequencies  (0.1-0. 5  GHz)  maximum  reflections  occur  at  the  top  and 
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Fig.  23.  Floe  size  distribution  obseirved  in  the  central  Beaufort 
Sea  during  mid-August  (based  on  measurements  by  Rothrock 
and  Thorndike ,  1984). 

bottom  of  the  ice,  providing  in  principle  a  means  to  obtain  thickness 
information.  At  higher  frequencies  (several  GHz)  radar  images  primar  ily 
show  differences  in  surface  roughness.  Passive  methods  measure  thermal 
emission  from  the  ice,  generally  in  the  3-100  GHz  range.  A  serious 
limitation  with  passive  observations  is  lack  of  spatial  resolution. 
Whereas  satellite -borne  radar  can  resolve  surface  elements  tens  of 


meters  across,  passive  resolution  varies  from  150  km  at  6  GHz  to  30  km 
at  37  GHz.  On  the  plus  side,  passive  techniques  can  distinguish  several 
different  ice  types  because  of  distinct  differences  in  thermal  emission. 
Understanding  microwave  signatures  from  various  types  of  ice  requires 
Information  on  the  dielectric  constant,  a  complex  quantity  whose  real 
part  describes  surface  reflection  and  backscattering  and  whose  Imaginary 
part  describes  attenuation  within  the  ice.  Dielectric  properties  of  sea 
ice  depend  on  the  properties  of  its  constituents  (ice,  brine,  air.  solid 
salts) ,  and  various  expressions  have  been  derived  to  express  bulk  values 
in  terms  of  temperature,  brine  volume,  and  pore  geometry  (I/eeks  and 
Ackley,  1982;  Morey  et  al.,  1984). 

The  availability  of  routine  microwave  data  from  polar  orbiting 
satellites  has  made  it  possible  to  procure  quantitative  Information  on 
ice  conditions  in  the  polar  oceans  (e.g. ,  Carsey,  1982;  Zwally  et  al., 
1983)  that  cannot  be  obtained  by  other  techniques.  Most  of  the  micro- 
wave  data  have  been  in  the  form  of  passive  "brightness  temperatures" 
which,  after  suitable  analysis,  can  give  ice  concentration,  ice  extent, 
and  the  relative  amounts  of  first-year  and  multiyear  ice.  Unfortu¬ 
nately,  the  procedure  for  translating  brightness  temperature  measure¬ 
ments  into  quantitative  estimates  of  ice  state  is  not  straightforward. 

A  single  brightness  temperature  is  a  spatial  Integral  of  energy  emitted 
from  a  region  that  is  large  (30-150  km  across)  compared  to  the  scale  of 
individual  features  such  as  leads  and  floes.  Even  though  such  features 
cannot  be  resolved  by  passive  microwave  imagery,  it  is  possible  to  esti¬ 
mate  their  contribution  to  an  observed  brightness  temperature  and  so 
infer  their  relative  abundance.  The  method  makes  use  of  frequency- 
dependent  differences  in  microwave  emissivity ,  the  ratio  of  energy 
actually  emitted  to  that  emitted  by  a  blackbody  under  identical  condi¬ 
tions.  In  the  Arctic,  multiyear  ice,  first-year  ice,  and  seawater 
exhibit  markedly  different  changes  in  emissivity  with  frequency  (Fig. 
24).  At  lower  frequencies  first-year  and  multiyear  ice  appear  much  the 
same,  while  leads  and  open  water  have  much  lower  brightness  temperatures 
owing  to  a  lower  emissivity.  At  higher  frequencies  first-year  ice 
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has  a  surface  layer  containing  little  brine,  but  many  air  pockets. 
Radiation  emitted  from  brine  deeper  in  the  ice  is  subject  to  scattering 
by  the  air  pockets  as  it  moves  upward  through  the  ice.  The  degree  of 
scattering  Increases  with  frequency,  and  hence  the  observed  decrease  in 
emissivity.  In  first-year  ice,  brine  is  present  very  close  to  the 
surface  so  that  the  emitted  radiation  undergoes  little  scattering,  and 
as  a  result  is  relatively  Independent  of  frequency.  It  is  thus  possible 
to  use  lower  frequency  data  to  determine  ice  concentration  and  extent, 
and  higher  frequency  data  to  calculate  first-year/multiyear  ratios.  The 
method  is  of  limited  value  during  the  summer  because  melting  erases  any 
distinction  between  ice  types  and  because  leads  cannot  be  distinguished 
from  surface  melt  ponds.  Ice  extent,  however,  can  still  be  found, 
regardless  of  season  or  the  location  of  the  edge. 

Uncertainties  inherent  in  the  determination  of  ice  concentration 
are  larger  than  the  amount  of  open  water  normally  present  within  the 
arctic  ice  pack  during  the  growing  season,  so  microwave  data  are 
primarily  used  to  map  the  distribution  of  first-year  and  multiyear  ice 
in  the  Northern  Hemisphere.  Because  of  weak  surface  melting  during  the 
summer,  differences  in  microwave  emission  from  first-year  and  multiyear 
ice  tend  to  be  subdued  or  totally  lacking  in  the  Antarctic,  particularly 
in  the  Weddell  Sea  region.  Although  signatures  similar  to  multiyear  ice 
have  been  seen  in  the  Ross,  Amundsen,  and  Bellingshausen  seas,  little  is 
yet  known  regarding  the  physical  characteristics  of  this  ice.  For  this 
reason,  microwave  analyses  of  ice  conditions  in  the  Southern  Hemisphere 
are  generally  focused  on  seasonal  and  spatial  variations  in  ice 
concentration  (e.g.,  Zwally  et  al.,  1983). 


5.  HEAT  AND  MASS  BALANCE 


Although  an  Ice  cover  reduces  momentvun  transfer  between  the  air  and 
the  water,  thermal  processes  associated  with  Its  growth  and  movement 
have  the  greatest  Impact  on  the  ocean.  In  particular,  solar  Input  to 
the  upper  ocean  Is  drastically  reduced,  limiting  biological  activity 
beneath  the  Ice.  The  cold  dense  brine  rejected  by  growing  Ice  promotes 
convective  overturning,  causing  any  heat  contained  In  the  mixed  layer  to 
be  transported  upward  where  It  Is  lost  to  the  Ice.  Thus,  unlike  the 
open  ocean,  heat  storage  In  the  upper  part  of  the  water  column  Is 
minimal.  Heat  and  mass  balance  measurements  have  been  carried  out  from 
land-based  stations  In  the  Arctic  and  Antarctic,  as  well  as  from  manned 
and  unmanned  drifting  stations  In  the  Arctic  Basin.  Conditions  In  the 
coastal  regions,  however,  are  Influenced  by  the  proximity  of  the  nearby 
land  masses  and  are  not  necessarily  representative  of  the  Ice  pack  as  a 
whole.  As  a  result,  our  understanding  of  Interactions  between  the  Ice, 
ocean,  and  atmosphere  Is  based  largely  on  what  has  been  seen  In  the 
perennially  Ice -covered  portions  of  the  Arctic  Basin.  The  extent  to 
which  such  Information  applies  to  more  dynamically  active  areas  such  as 
the  Marginal  Ice  Zone  and  the  antarctic  Ice  pack  Is  unknown. 

A.  Energy  Fluxes 
I .  Radiation 

Radiation  dominates  heat  exchange  In  the  polar  regions  throughout 
the  year.  All  matter  radiates  energy  in  proportion  to  the  fourth  power 
of  its  absolute  temperature,  i.e.,  F  -  eoT^  ,  where  F  is  the  flux  of 
emitted  radiation,  c  Is  the  surface  emlsslvlty,  and  o  Is  the  Stefan- 
Boltzmann  constant  -  5.67'10*®  W  m'^  °K’^(4.88'10’^  cal  cm*^  s'^  °K‘^) . 
Most  of  the  radiation  reaching  the  earth's  surface  lies  in  two  spectral 
bands:  shortwave  radiation ,  which  originates  from  the  sun  and  has  a 
narrow  intensity  maximum  near  X  -  500  nm,  and  longwave  radiation,  which 
originates  from  the  earth  and  atmosphere  with  a  broad  Intensity  maximum 
near  X  -  12000  nm.  Because  there  is  little  overlap  in  the  spectral 
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distribution  of  radiation  derived  from  these  two  sources,  there  Is  a 
natural  separation  In  their  measurement  and  treatment. 

Downward  longwave  radiation  (Fj^)  emitted  from  clouds  and  water 
vapor  In  the  atmosphere  Is  a  major  component  In  the  surface  heat 
balance.  During  most  of  the  winter  and  fall  Incoming  shortwave  radi¬ 
ation  (F^)  Is  negligible  and,  even  at  the  height  of  the  melt  season,  F^. 
Is  only  slightly  larger  than  F^  (Table  1).  On  an  annual  basis  F^  Is 
more  than  double  F^  In  the  Central  Arctic.  Despite  the  magnitude  of  Fj^, 
the  net  effect  of  longwave  radiation  Is  to  remove  heat  from  the  Ice. 

This  Is  because  the  radiative  temperature  of  the  Ice  usually  exceeds 
that  of  the  atmosphere  so  that  longwave  radiation  emitted  by  the  Ice 
(F'*')  Is  almost  always  greater  than  Fj^.  Measurements  of  F^  and  F'*"  are 
difficult  In  the  polar  regions  because  of  constant  problems  with  dew  and 
frost  formation  on  the  Instruments.  Fortunately,  both  can  be  estimated 
from  other  Information.  Surface  temperature  data  are  frequently  used  to 
calculate  F^  according  to  the  formula,  F^  -  where  Is  the 

longwave  emlsslvlty  which  equals  0.99  for  snow  and  0.97  for  bare  Ice  and 
melt  ponds.  Because  a  large  part  of  the  Incoming  longwave  radiation  Is 
emitted  by  clouds  and  water  vapor  close  to  the  ground,  F^  can  be 
approximated  fairly  well  from  surface  observations  of  air  temperature 


Table  1.  Monthly  and  annual  radiation  totals  (MJ  m"^)  over  perennial  ice  in  the 
Central  Arctic  (after  Marshunova,  1961;  Maykut  and  Untersteiner,  1971) 
Numbers  in  parentheses  Include  the  effects  of  summer  melt  ponds . 


JAN 

FEB 

MAR 

APR 

MAY 

JUN 

JUL 

AUG 

SEP 

OCT 

NOV 

DEC 

YEAR 

Incoming  longwave 
radiation  [F^] 

435 

430 

430 

495 

630 

750 

800 

780 

690 

580 

475 

450 

6945 

Outgoing  longwave 
radiation  [F^] 

-520 

-500 

-505 

-600 

-725 

-830 

-835 

-820 

-720 

-610 

-535 

-520 

-7720 

Net  longwave 
radiation  (Fj^+F^l 

-85 

-70 

-75 

-105 

-95 

-80 

-35 

-40 

-30 

-30 

-60 

-70 

-775 

T  coming  shortwave 
radiation  (F^l 

100 

420 

730 

785 

575 

370 

155 

20 

3155 

Absorbed  shortwave 
radiation  ((l-o)F^) 

15 

80 

145 

170 

245 

(295) 

115 

(165) 

35 

5 

810 

(910) 

Net  radiation 

-85 

-70 

-60 

-25 

50 

90 

210 

(255) 

75 

(120) 

5 

-25 

-60 

-70 

35 

(125) 

(T^) ,  partial  pressure  of  water  vapor  (e^) ,  and  fractional  cloud  cover 
(C) .  Longwave  paraneterlzatlons  have  the  form,  Fj^  -  e*aT^^,  where  e*  is 
the  effective  emlssivlty  of  the  atmosphere.  Several  empirical 
expressions  have  been  suggested  for  e*  in  the  Arctic.  The  simplest, 

e*  -  0.7829(1  +  0.2232c2-^^)  (24) 

obtained  by  Maykut  and  Church  (1973)  from  observations  at  Point  Barrow, 
Alaska,  gives  estimates  of  accurate  to  within  about  7%  under  overcast 
skies  and  12%  under  clear  skies.  No  correlation  was  found  between  e 
and  e^^.  More  complex  expressions  depending  on  e^,  C,  and  a 
seasonally  varying  empirical  parameter  have  been  obtained  by  Marshunova 
(1961)  and  Idso  and  Jackson  (1969). 

The  sun  being  continuously  above  or  below  the  horizon  for  periods 
of  2-6  months  at  a  time  produces  seasonal  changes  in  F^.  that  are  much 
larger  than  those  in  F^.  In  spite  of  extensive  clouds  and  low  sun 
angles,  monthly  totals  of  F^.  in  the  Arctic  can  exceed  those  measured  at 
many  midlatitude  stations  (e.g..  New  York,  Chicago)  during  the  summer. 
Spatial  variations  in  F^  poleward  of  65°N  tend  to  be  small  during  the 
summer  months  (Marshunova,  1961).  Because  of  the  extreme  summer 
cloudiness  (C  -  0.8-0. 9),  most  of  the  shortwave  radiation  reaching  the 
surface  is  diffuse  rather  than  direct.  Clouds  over  the  polar  oceans 
tend  to  be  thin  and  so  typically  reduce  F^  by  30-50%  as  compared  to  the 
80-90%  often  observed  at  lower  latitudes.  The  relatively  low  water 
vapor  content  of  arctic  clouds  allows  substantial  amounts  of  near 
Infrared  energy  to  reach  the  surface.  Under  clear  skies  about  55%  of  F^ 
lies  outside  the  visible  region  (400-700  nm) ,  while  under  overcast  skies 
this  figure  decreases  to  about  35%.  Shortwave  parameterizations  used  in 
large-scale  heat  balance  studies  and  numerical  models  are  generally 
written  in  the  form  F^  -  a  where  F^^  is  the  incoming  shortwave 

radiation  at  the  surface  in  the  absence  of  clouds  and  a  is  the 
fractional  attenuation  by  clouds.  Although  both  F^^  and  a  can  be 
calculated  theoretically,  the  procedures  are  too  involved  for  many 


applications  and  empirical  relationships  are  frequently  employed.  A 
particularly  convenient  form  for  was  given  by  Zillman  (1972), 
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a 


(25) 


where  is  the  solar  constant,  e^  is  in  millibars  and  z  is  the  solar 
zenith  angle,  a  function  of  latitude,  day,  and  hour.  Limited  tests  with 
arctic  data  (Maykut,  in  press)  suggest  that  (25)  tends  to  underestimate 
at  high  latitudes.  The  reason  is  probably  related  to  the  smaller 
amount  of  water  vapor  in  the  polar  atmosphere  and  to  backscattering  of 
radiation  reflected  from  the  ice.  This  suggests  that  a  high  latitude 
correction  may  be  desirable  for  formulas  based  on  observations  at  lower 
latitudes.  On  the  basis  of  model  calculations.  Shine  (1984)  reports 
that  improved  performance  can  be  obtained  from  (25)  over  ice  and  snow  if 
the  denominator  is  changed  to  1.2  cos  z  +  (1.0  +  cos  z)ej^’10  +  0.046. 

Expressions  for  a  have  been  obtained  by  Marshimova  (1961)  [a  -  1  -  kC, 
where  k  increases  from  0.15  in  March  to  0.5  in  July]  and  Laevastu  (1960) 
[a  -  1  -  0.6C  ].  A  parameterization  for  the  net  shortwave  flux  over 
ice  and  snow,  claimed  to  be  accurate  to  within  a  few  percent  under 
cloudy  arctic  conditions,  has  also  been  derived  by  Shine  (1984)  in  terms 
of  solar  zenith  angle,  surface  albedo,  and  the  optical  thickness  of  the 
cloud  cover. 

As  mentioned  earlier,  the  summer  melt  cycle  over  perennial  ice  is 
driven  almost  entirely  by  absorbed  shortwave  radiation.  Although  F^. 
reaches  a  maximum  in  June,  the  largest  values  of  net  shortwave  radiation 
occur  in  July  because  of  the  decrease  in  a  following  removal  of  the  snow 
cover.  The  presence  of  melt  ponds  significantly  increases  the  net 
shortwave  radiation  and  hastens  the  rate  of  decay  and  retreat  of 
seasonal  ice.  In  areas  of  perennial  ice,  ponds  increase  shortwave  input 
to  the  ice  by  20-40%  (Table  1),  but  have  only  a  small  effect  on  the 
annual  average  ice  thickness  (Maykut  and  Untersteiner,  1971).  Ice 
concentration  is  also  an  Important  parameter  because  leads  admit  large 
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quantities  of  energy  to  the  ocean.  Once  the  snow  cover  has  vanished, 
the  amount  of  energy  available  for  primary  production  in  the  underlying 
water  Is  largely  determined  by  the  relative  area  occupied  by  leads  and 
melt  poiids .  Table  2  compares  how  1  m  thick  layers  of  melt  pond/blue  ice 
and  white  ice  distribute  for  different  ice  concentrations.  When 
is  high,  transmission  through  blue  ice  accounts  for  most  of  the  energy 
absorbed  beneath  the  ice .  As  A^  decreases ,  leads  quickly  begin  to 
dominate  solar  input  to  the  water  beneath  the  ice. 


Table  2.  The  effect  of  ice  concentration  and  surface  type  on  absorp¬ 
tion  and  transmission  of  incoming  shortwave  radiation  in  a 
region  of  1  m  thick  bare  ice.  Values  are  expressed  as  per¬ 
centages  of  weighted  by  relative  area. 


Ice 

Surface 

Net  input 

Z  Absorbed 

Z  Transmitted 

Regional 

concentration 

type 

at  surface 

0  £  z  £  1  m 

below  1  m 

input 

blue  ice 
white  ice 


leads 


blue  ice 
white  ice 


leads 


2.  Turbulent  Exchange 

Apart  from  outgoing  longwave  radiation,  the  only  means  of 
transferring  energy  from  the  ice  to  the  atmosphere  is  through  turbulent 
exchange.  The  sensible  heat  flux  (Fg)  results  from  temperature 
differences  between  the  ice  and  atmosphere,  while  the  latent  heat  flux 
(Fg)  is  produced  by  the  release  of  latent  heat  associated  with 
evaporation,  condensation,  and  sublimation.  These  fluxes  are  frequently 
dismissed  as  being  unimportant  because  they  are  typically  1-2  orders  of 
magnitude  smaller  than  the  radiation  fluxes.  This  argument  is 
misleading,  however,  because  the  net  radiation  balance  is  small  in 
comparison  to  its  individual  components  (Table  1).  In  fact,  the 
magnitude  of  Fg  +  Fg  in  the  Central  Arctic  is  on  the  order  of  20-50%  of 


the  total  net  radiation.  Theoretical  calculations  (Haykut  and 
Unter Steiner ,  1969)  have  demonstrated  that  the  equilibrium  thickness  of 
multiyear  ice  is  very  sensitive  to  changes  in  Fg  and  Fg  during  the 
svunmer.  Turbulent  fluxes  near  the  boundaries  of  the  ice  pack  tend  to  be 
larger  and  appear  to  be  a  significant  factor  in  the  summer  decay  cycle 
there . 

The  rate  of  turbulent  heat  transfer  depends  on  surface  roughness , 
wind  speed,  boundary  layer  stability,  and  the  gradients  in  temperature 
and  water  vapor  near  the  surface.  Commonly  used  parameterizations  for 
these  fluxes  are:  Fg  -  PaCpCs«(Ta  *  T^)  and  Fg  -  Pa^vCgUCqg  -  q^) , 
where  is  the  density  of  the  air,  Cp  is  the  specific  heat  of  the  air, 
u  is  the  wind  speed,  is  the  latent  heat  of  vaporization,  q^  and  q^ 

are  specific  humidities  in  the  air  and  at  the  surface,  and  C.  and  C_  are 
bulk  transfer  coefficients  which  vary  with  stability  and  surface  rough¬ 
ness  .  Eddy  correlation  measurements  over  multiyear  ice  (Thorpe  et  al. , 
1973)  indicate  that  Cg  -  1.2' 10*^  and  Cg  -  0.55*10’^,  while  observations 
in  the  unstable  boundary  layer  over  a  refreezing  lead  (Lindsay,  1976) 
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give  values  of  2. 0-2. 5 '10  for  both  Cg  and  Cg.  Differences  in  surface 
geometry  (ridges,  hummocks,  floe  sizes)  can  be  expected  to  produce  local 
variations  in  Cg  and  Cg  of  at  least  a  factor  of  2. 

Studies  of  turbulent  heat  exchange  over  perennial  sea  ice  (Doronin, 
1963;  Badgley,  1966;  Leavitt  et  al.,  1978)  show  that,  unlike  the 
situation  over  temperate  and  tropical  oceans,  Fg  is  negligible  during 
much  of  the  polar  year  (November-April  in  the  Arctic) .  Even  during  the 
summer,  latent  heat  losses  are  only  slightly  larger  than  those  from 
sensible  heat.  During  the  winter  when  there  is  a  persistent  temperature 
inversion  close  to  the  surface,  Fg  supplies  heat  to  the  ice,  compensat¬ 
ing  for  one- third  to  one-half  the  net  longwave  loss.  Although  daily 
totals  of  Fg  during  the  summer  are  comparable  in  magnitude  to  the  winter 
values,  monthly  averages  are  small  as  a  result  of  frequent  alternations 
in  the  direction  of  Fg.  All  studies  in  the  Central  Arctic  show  a  net 
turbulent  heat  loss  for  the  summer  that  acts  to  reduce  the  amount  of 
surface  melting.  Without  these  turbulent  heat  losses,  surface  ablation 


would  be  20-100%  larger  and  the  equilibrium  thickness  reduced  by  up  to 
50%  (Maykut  and  UntersCeiner,  1969).  At  the  margins  of  the  ice  pack,  Fg 
and  Fg  are  generally  downward  and  contribute  substantial  amounts  of  heat 
to  the  surface.  During  May  and  June,  for  example,  the  turbulent  fluxes 
can  account  for  as  much  as  one -third  of  the  net  heat  input  to  coastal 
ice  (Langleben,  1966). 

3.  Heat  Conduction 

Even  during  the  melt  season,  there  is  an  upward  conduction  of  heat 

in  the  ice  [F^.  -  kj^(3T/32)]  which  supplies  much  of  the  energy  needed  to 

balance  the  longwave  loss  at  the  surface.  In  thin  ice  F-(z  -  0)  =  F 

c  c 

(z  -  h) ,  greatly  simplifying  the  treatment  of  young  ice  growth.  In 
thicker  ice  Internal  phase  changes  and  heat  storage  become  important 
and  growth  calculations  must  take  into  account  nonlinearities  in  the 
temperature  profile.  Over  an  annual  cycle  the  amount  of  heat  conducted 
to  the  surface  of  multiyear  ice  is  the  sum  of  the  sensible  heat  input 
from  the  ocean  (F^^) ,  the  latent  heat  released  by  net  bottom  accretion, 
and  the  shortwave  energy  temporarily  stored  in  brine  pockets  and  surface 
melt  ponds.  For  perennial  arctic  ice  F^.  *  6.6  kcal  cm*^  year'^ 

(8.8  W  m*^),  of  which  23%  is  supplied  by  the  ocean,  39%  by  bottom 
accretion,  29%  by  release  of  stored  heat  in  the  brine  pockets,  and  9%  by 
the  refreezing  of  surface  meltwater  (Maykut,  in  press).  F^.  can  be  more 
than  two  orders  of  magnitude  larger  in  young  ice  than  in  multiyear  ice. 
Over  the  course  of  the  entire  growth  season  the  amount  of  heat  conducted 
to  the  surface  of  seasonal  ice  is  roughly  2-3  times  the  annual  total  for 
perennial  ice. 

4.  Oceanic  Heat  Flux 

A  variety  of  observational  and  theoretical  evidence  points  to  a 
significant  input  of  heat  from  the  ocean  to  the  ice,  but  direct  measure¬ 
ments  of  this  flux  (F^^)  have  not  been  successful.  In  the  Antarctic 
where  surface  ablation  is  small,  the  seasonal  disappearance  of  the  ice 
must  be  caused  almost  entirely  from  energy  contained  in  the  water. 


Certainly  shortwave  radiation  absorbed  In  the  upper  ocean  Is  a  major 
component  of  there;  however,  heat  from  the  deeper  ocean  is  also 
important.  Gordon  (1981)  estimates  that  up  to  50%  of  could  be 
derived  from  water  below  the  pycnocline,  but  the  uncertainty  in  this 
estimate  appears  to  be  large.  In  the  Central  Arctic  F„  is  about 
1.5  kcal  cm*  year*  (2  W  m'^).  The  source  of  this  heat  was  originally 
believed  to  be  the  layer  of  relatively  warm  Atlantic  water  that  occurs 
across  much  of  the  Arctic  Basin.  More  recently  it  has  been  suggested 
that  temperature  changes  in  the  Atlantic  layer  are  largely  due  to  mixing 
with  colder  shelf  water  (Aagaard  et  al.,  1981)  and  that  F^  is  determined 
almost  entirely  by  shortwave  input  through  leads  and  areas  of  thin  ice 
(Maykut,  1982).  In  areas  of  seasonal  ice  the  upper  part  of  the  water 
column  at  the  onset  of  freezing  can  contain  1-2  kcal  cm*  (40-80  MJ  m*  ) 
which  is  slowly  lost  as  the  mixed  layer  deepens  in  response  to  salt 
rejection  by  the  growing  ice.  The  largest  values  of  F^  in  the  Arctic 
probably  occur  in  the  Fram  Strait  region  where  the  Atlantic  enters  the 
basin  and  is  still  close  to  the  surface. 


B.  Response  of  the  Ice  to  Environmental  Changes 

Because  sea  ice  appears  to  be  an  integral  part  of  the  global 
climate  system,  there  are  numerous  questions  that  need  to  be  answered 
regarding  its  long-term  response  to  various  environmental  factors.  For 
example,  how  will  the  polar  ice  pack  be  affected  by  the  increasing 
concentration  of  carbon  dioxide  in  the  atmosphere?  Can  man  artificially 
alter  ice  conditions  to  limit  these  effects?  How  does  the  state  of  the 
ice  pack  change  over  a  glacial  cycle?  Although  some  information  on 
historical  changes  in  the  ice  pack  can  be  gained  from  the  fossil  record 
and  oxygen  Isotope  analysis  of  ocean  sediments,  sensitivity  studies  with 
coupled  ocean- ice -atmosphere  models  are  needed  to  fully  understand  the 
response  of  the  ice.  Such  models  are,  however,  still  in  a  relatively 
crude  state  of  development  and  do  not  adequately  simulate  many  system 
interactions  and  feedbacks  that  are  likely  to  affect  ice  extent  and 
thickness . 


Even  though  prediction  of  future  changes  in  the  state  of  the  ice 
cannot  yet  be  made  with  any  real  degree  of  confidence,  it  is  possible  to 
determine  with  some  accuracy  how  specific  changes  in  such  factors  as  air 
temperature,  snowfall,  cloudiness,  and  oceanic  heat  flux  will  alter  ice 
thickness.  This  is  the  first  step  needed  in  addressing  the  more 
difficult  problem  of  ice  extent.  To  estimate  thickness  changes,  it  is 
first  necessary  to  formulate  energy  balance  equations  at  the  upper  and 
lower  boundaries  of  the  ice.  The  flux  balance  at  the  surface  is 

BT^ 

(1  -  ct)(l  -  +  Fl  -  F+  +  F3  +  Fe  +  =  0  .  (26) 


where  Fjjj  -  p j^L^Cdh/dt)^  -  the  heat  loss  due  to  surface  melting.  When  T^ 
is  below  freezing,  Fj^^  -  0;  when  T^  is  at  the  melting  point,  any  surplus 
of  energy  will  be  balanced  by  melting  and  a  change  in  h.  At  the  bottom 
of  the  ice 


Hh  ^'^’l 


(27) 


Equations  (26)  and  (27)  are  coupled  via  the  heat  conduction  terms, 
ki(3Ti/9z)o  and  kj^(3T£/3z)j^,  whose  calculation  requires  knowledge  of  the 
vertical  temperature  distribution  within  the  ice.  Changes  in  T^(z)  can 
be  calculated  using  a  form  of  the  thermal  diffusion  equation  modified  to 
Include  the  effects  of  internal  heating  from  solar  radiation 


3Ti 


3T^ 


11  3t 


-Ki(z-Zo) 


3z 


(28) 


Simultaneous  solution  of  (26) -(28)  then  makes  it  possible  to  determine 
the  thermal  response  of  the  ice  to  any  of  the  parameters  explicitly  or 


implicitly  contained  in  the  equations.  Inclusion  of  a  snow  cover  adds 
an  additional  boundary  condition  and  heat  conduction  equation,  but  does 
not  significantly  complicate  the  solution.  Large-scale  simulations  of 
the  ice  cover  generally  retain  (26)  and  (27),  but  simplify  the  problem 
either  by  eliminating  (28)  through  the  assumption  of  a  linear 
temperature  profile  in  the  ice,  or  by  resolving  Tj^(z)  only  crudely. 

Model  results  (Haykut  and  Untersteiner ,  1971)  indicate  that  the 
equilibrium  thickness  (see  Section  2B)  is  especially  sensitive  to 
factors  that  affect  the  amount  of  summer  melting:  air  temperature, 
surface  albedo,  incoming  radiation,  and  turbulent  heat  exchange. 

Complete  disappearance  of  perennial  arctic  ice  would,  for  example,  occur 
with  increases  of  3-5°C  in  T^,  25-30%  in  incident  solar  radiation,  or  a 
15-20%  decrease  in  summer  albedo.  Changes  in  cloudiness  affect  both  the 
shortwave  and  longwave  radiation.  Decreased  cloudiness  lowers  the 
amount  of  incoming  longwave  radiation  while  increasing  the  downward 
shortwave  flux.  A  decrease  in  cloudiness  thus  acts  to  cool  the  ice 
during  the  winter,  but  has  only  a  weak  effect  during  the  summer  because 
the  opposing  changes  in  incident  radiation  tend  to  balance.  Year  to 
year  variations  in  arctic  cloudiness  are  presently  on  the  order  of  25% 
and  do  not  appear  to  change  hg  significantly. 

Theoretical  predictions  of  how  variations  in  snowfall  and  oceanic 
heat  flux  affect  hg  in  the  Central  Arctic  are  shown  in  Figs.  25  and  26. 
In  spite  of  the  large  effect  that  snow  depth  has  on  ice  temperatures  and 
growth  rates,  differences  in  snowfall  are  found  to  have  little  impact  on 
hg  until  depths  exceed  about  twice  the  present  day  average  of  40  cm. 

With  less  snow  the  ice  is  colder  and  winter  growth  greater,  but  summer 
ice  melting  is  also  greater  because  of  earlier  removal  of  the  snow. 
Changes  in  surface  melting  and  bottom  accretion  approximately  balance 
until  snow  depths  become  so  large  that  ice  ablation  is  delayed  until 
late  in  the  summer  when  begins  to  decrease  rapidly  with  time.  Above 
hg  -  80  cm  the  decrease  in  ice  ablation  is  greater  than  the  decrease  in 
accretion,  and  hg  increases  steeply  with  increasing  hg .  Under  present 
surface  conditions  the  thickness  of  undeformed  perennial  ice  will  not 
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Fig.  25.  Average  equlllbrluo  thickness  of  perennial  arctic  sea 
Ice  as  a  function  of  Daxlmum  annual  snow  depth  (after 
Haykut  and  Untersteiner ,  1971) . 


exceed  7  m  In  the  Central  Arctic,  even  In  regions  where  the  oceanic  heat 
flux  Is  negligible  (Fig.  26).  However,  If  h^  Is  In  excess  of  about  120 
cm,  summer  melting  In  the  Arctic  Is  generally  Insufficient  to  remove  the 
snow  cover  completely,  and  equilibrium  Ice  will  grow  from  above  and 
ablate  from  below,  becoming  In  effect  a  "floating  glacier."  The 
combination  of  small  F^  and  high  snowfall  can  produce  Ice  of  almost  any 
thickness  and  has  been  used  to  explain  the  origin  of  very  thick  Ice 
found  In  Nansen  Sound  and  along  the  north  coast  of  Ellesmere  Island 


EQUILIBRIUM  ICE  THICKNESS  (m) 


Fig.  26.  Equilibrium  thickness  of  arctic  sea  Ice  as  a  function  of 
the  average  annual  oceanic  heat  flux.  Average  annual 
thickness,  as  well  as  the  absolute  annual  maximum  and 
minimum,  are  shown  (after  Maykut  and  Untersteiner, 

1971). 

(Walker  and  Wadhams,  1979).  The  calculations  also  predict  that  the 

arctic  Ice  cover  would  be  seasonal  If  totaled  more  than  about  5.5 
-2  -1  -2 

kcal  cm  year  (7.3  W  m  ).  It  is  suspected  that  average  values  of  F^ 
In  many  parts  of  the  Southern  Ocean  may  exceed  this  value  by  at  least  a 
factor  of  2-3  (Gordon  et  al.,  1984). 


It  was  shown  In  Fig.  5  that,  for  constant  thermal  forcing,  winter 
growth  rates  In  young  Ice  are  strongly  dependent  on  thickness.  The  same 
holds  true  for  the  surface  heat  balance.  Surface  temperatures  In  thin 
Ice  decrease  rapidly  with  Increasing  thickness  (Fig.  27),  affecting  the 
conductive  heat  flux,  the  turbulent  exchange,  and  the  longwave  loss.  In 


- 


shortwave  absorption  as  the  ice  grows.  Heat  exchange  over  thin  Ice 
tends  to  be  much  more  vigorous  than  over  thicker  Ice.  Turbulent  heat 
losses  over  a  refreezlng  lead,  for  example,  can  be  more  than  two  orders 
of  magnitude  larger  than  those  over  nearby  thick  Ice.  Thus,  relatively 
small  areas  of  thin  Ice  and  open  water  can  have  a  significant  Impact  on 
regional  heat,  mass,  and  salt  balances. 
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Theoretical  calculations  have  been  carried  out  to  quantify  how 
thickness  affects  the  surface  heat  balance  In  regions  of  perennial 
arctic  Ice  {Haykut,  1978).  Results  are  svunmarlzed  In  Table  3.  When  the 
Ice  Is  thin,  heat  Input  at  the  surface  Is  controlled  by  F^,  and 
temperatures  are  relatively  Insensitive  to  changes  In  atmospheric 


Table  3.  Surface  heat  balance  components  In  the  Central  Arctic  as  a  function  of 
ice  thickness  (after  Mayloit,  1978).  Energy  fluxes  are  given  In  W  m~^ 
and  apply  to  conditions  on  the  first  day  of  the  month. 


Sept. 

Oct. 

Nov.  Dec. 

Jan. 

Feb. 

March 

April 

May 

June 

Open  H'aitr 

Net  thonwave  radiation 

S9 

24 

7 

83 

209 

281 

Net  longwave  radiation 

-21 

-69 

-114 

-138 

-142 

-147 

-149 

-140 

-98 

-43 

Sensible  heat  flux 

-68 

-259 

-4S8 

-540 

-543 

-575 

-615 

-520 

-276 

-33 

Latent  heat  flux 

-34 

-108 

-139 

-145 

-145 

-147 

-ISO 

-144 

-112 

-31 

Oceanic  heat  flux 

130 

436 

711 

823 

830 

869 

914 

804 

486 

107 

OOSm 

Available  shortwave  energy* 

60 

16 

5 

56 

141 

189 

Net  longwave  radiation 

-22 

-42 

-70 

-87 

-91 

-93 

-93 

-93 

-76 

-51 

Sensible  heat  flux 

-44 

-142 

-261 

-311 

-312 

-334 

-362 

-313 

-179 

-74 

Latent  heat  flux 

-31 

-so 

-49 

-47 

-46 

-45 

-44 

-51 

-62 

-47 

Conductive  heat  flux 

37 

218 

380 

445 

449 

472 

494 

401 

176 

-17 

0.1  m 

Available  shortw  ave  energy 

S6 

IS 

4 

52 

131 

ITS 

Net  longwave  radiation 

-20 

-33 

-55 

-69 

-73 

-75 

-73 

-77 

-68 

-50 

Sensible  heat  flux 

-36 

-101 

-186 

-222 

-222 

-238 

-262 

-232 

-143 

-73 

Latent  heat  flux 

-27 

-34 

-30 

-27 

-26 

-25 

-24 

-31 

-47 

-47 

Conductive  heat  flux 

27 

IS3 

271 

318 

321 

338 

355 

288 

127 

-5 

0.2  m 

Available  shortwave  energy 

S3 

14 

4 

49 

124 

166 

Net  longwave  radiation 

-19 

-25 

-41 

-53 

-57 

-58 

-55 

-62 

-61 

-49 

Sensible  heat  flux 

-29 

-64 

-115 

-137 

-136 

-147 

-165 

-154 

-109 

-68 

Latent  heat  flux 

-23 

-22 

-16 

-14 

-14 

-13 

-12 

-17 

-34 

-46 

Conductive  hett  flux 

18 

97 

172 

204 

207 

218 

228 

184 

80 

-3 

0.4  m 

Available  shortwave  energy 

48 

13 

4 

46 

114 

153 

Net  longwave  radiation 

-17 

-19 

-31 

-41 

-45 

-45 

-42 

-51 

-55 

-47 

Sensible  heat  flux 

-22 

-36 

-62 

-72 

-70 

-77 

-91 

-93 

-82 

-60 

Latent  heat  flux 

-19 

-14 

-9 

-6 

-6 

-6 

-5 

-10 

-26 

-44 

Conductive  heat  flux 

10 

S6 

102 

119 

121 

128 

134 

108 

49 

-2 

0.8  m 

A  vailable  shortwave  energy 

4S 

12 

3 

42 

104 

140 

Net  longwave  radiation 

-16 

-16 

-24 

-33 

-37 

-37 

-34 

-43 

-51 

-46 

Sensible  heat  flux 

-17 

-18 

-27 

-» 

-27 

-31 

-42 

-53 

-61 

-54 

Latent  heat  flux 

-16 

-9 

-4 

-3 

-2 

-2 

-2 

-5 

-19 

-39 

Conductive  heat  flux 

4 

31 

55 

66 

66 

70 

75 

59 

27 

-1 

3.0  m 

Available  shortw  ave  energy 

16 

4 

1 

17 

42 

59 

Net  longwave  radiation 

-14 

-12 

-18 

-26 

-31 

-29 

-25 

-32 

-38 

-35 

Sensible  heat  flux 

-S 

0 

6 

II 

17 

16 

12 

9 

-2 

-8 

Latent  heat  flux 

-9 

-4 

0 

0 

0 

0 

0 

0 

-4 

-10 

Conductive  heat  flux 

12 

12 

12 

IS 

14 

13 

12 

6 

2 

-6 

forcing.  Heat  exchange  with  the  atmosphere  over  young  ice  is  dominated 
by  the  sensible  heat  flux.  In  March,  for  example,  ranges  from  a 
small  heat  gain  over  thick  ice  to  a  loss  over  open  leads  that  is  about 
3.5  times  the  total  of  the  incident  radiation  flxixes.  While  over 
open  leads  is  comparable  to  the  radiation  fluxes  during  most  of  the 
winter,  its  Importance  decreases  rapidly  with  h  owing  to  a  steep  decline 
in  saturation  vapor  pressure  below  0°C.  Over  thicker  first-year  ice  F^ 
is  comparable  to  the  net  longwave  radiation.  Energy  needed  to  support 
the  large  energy  fluxes  over  young  ice  is  derived  largely  from  latent 
heat  released  by  the  growth  of  ice,  and  not  by  heat  transport  in  the 
ocean.  In  all  cases  the  relative  differences  in  heat  exchange  tend 
toward  zero  in  the  spring  as  the  ice  approaches  the  melting  point.  The 
calculations  also  demonstrate  that  the  surface  heat  balance  becomes 
insensitive  to  thickness  once  h  exceeds  80-100  cm. 

D.  KeetbMl  Flaxes 

Large  changes  in  the  amount  of  thin  ice  and  open  water  can  occur  on 
time  scales  of  hours  to  days,  producing  corresponding  changes  in  the 
regional  heat  and  mass  balance  which  are  independent  of  the  incident 
energy  fluxes.  For  obvious  reasons  most  field  data  have  been  gathered 
from  thick,  relatively  stable  sea  ice.  However,  because  they  do  not 
take  into  account  contributions  made  by  areas  of  thin  ice  and  open 
water,  such  data  are  not  necessarily  representative  of  the  region  as  a 
whole.  The  large  spatial  variability  in  local  values  makes  direct 
determination  of  regional  fluxes  difficult.  An  alternate  approach  is  to 
utilize  theoretical  models  to  infer  large-scale  values  from  routine 
local  observations.  The  method  works  as  follows.  Let  be  any  property 
that  depends  on  h  (e.g.,  the  salt  flux  to  the  ocean).  The  contribution 
made  by  a  particular  thickness  category  is  found  by  multiplying  the 
fractional  area  covered  by  that  category  times  the  magnitude  of  ^  over 
that  thickness.  Summing  over  all  possible  thicknesses  then  gives  a 
total  for  the  region,  i.e.. 


♦(h)g(h)dh 


<p  is  defined  to  be  the  regional  or  large-scale  average.  In  practice  It 
Is  usually  necessary  to  work  with  discrete  thickness  categories.  If 
g(h)  Is  partitioned  Into  n  categories,  (29)  can  be  approximated  by 


7  ^  ♦iAGi  , 


where  <l>^  is  the  average  value  of  4i  In  the  category,  AG^  -  gj^Ah^, 


Ahf  “  hj^  -  h£_^,  and  g^  Is  the  average  value  of  g  In  the  1^"  category. 


Since  (p  and  g  are  also  usually  functions  of  time,  we  use  the  notation 
<4'>  -  Jo4'(t)dt  to  denote  spatial  averages  Integrated  over  time. 


Application  of  (29)  or  (30)  requires  general  heat  balance  and  Ice 
movement  Information.  A  thermodynamic  Ice  model  Is  used  with  heat 
balance  data  to  obtain  growth  rates  and  energy  fluxes  over  all 
thicknesses  of  Ice.  Growth  rates  and  Ice  velocities  are  then  used  to 
drive  an  Ice  thickness  distribution  model  and  calculate  g(h)  for  the 
region,  thereby  allowing  the  determination  of  from  (30). 

A  study  of  the  effects  of  thickness  variations  on  regional  fluxes 
has  been  carried  out  using  data  from  the  central  Beaufort  Sea  (Maykut, 
1982).  Results  are  summarized  In  Table  4  and  Figs.  28  and  29.  The 
calculations  clearly  demonstrate  that  relatively  small  amounts  of  thin 
Ice  can  play  a  major  role  In  the  overall  heat  and  mass  balance  of  a 
region.  Between  October  and  April  refreezlng  leads  generally  covered  no 
more  than  0.5%  of  the  area,  while  thin  Ice  (10  <  h  <  40  cm)  covered  only 
2*3%.  Greater  divergence  and  melting  of  thin  Ice  caused  a  large 
Increase  In  the  amount  of  open  water  during  the  spring  and  summer. 
Freezing  of  the  open  water  In  late  August  produced  a  September  peak  In 
the  distribution  of  young  Ice,  but  ridging  and  continued  growth  combined 
to  move  most  of  this  Ice  Into  the  80'*'  cm  category  by  October.  The 
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Fig.  28.  Calculated  distribution  of  ice  thickness  in  the  central 
Beaufort  Sea  (after  Maykut,  1982). 

annual  net  radiation  balance  for  the  region  was  more  than  double  what 
would  have  been  expected  on  the  basis  of  measurements  over  thicker 
(h  >  1  m)  ice.  Most  of  this  additional  energy  input  came  from  solar 
radiation  absorbed  in  summer  leads  and  ultimately  went  to  lateral 
melting  on  floe  edges  and  to  F^.  The  amount  of  solar  radiation  absorbed 
beneath  the  ice  during  the  summer  was  about  105  HJ  m*  ,  more  than  a 
factor  of  50  greater  than  the  amount  that  would  be  absorbed  beneath  a 
continuous  3  m  thick  ice  cover.  Conduction  of  heat  to  the  surface 
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Fig.  29.  Regional  Ice  production  In  the  central  Beaufort  Sea: 

(a)  annual  totals  In  each  thickness  category,  (b) 
monthly  totals  taking  Into  account  Ice  thickness  varla 
tlons  (stippled)  and  assuming  a  uniform  3  m  thick  Ice 
cover  (cross-hatched)  (after  Haykut,  1982). 


♦  I*! 


Increased  by  a  factor  of  2  when  thickness  variations  were  taken  Into 

account.  Most  of  this  heat  was  lost  to  the  atmosphere  during  the  fall 

and  winter  via  larger  turbulent  fluxes.  Turbulent  heat  Input  to  the 

_2 

atmosphere  totaled  about  200  HJ  m  for  the  year,  due  In  large  part  to 
greater  sensible  heat  losses  associated  with  young  Ice  formation  In  the 
fall.  It  Is  notable  that  the  open  water  category  (0-10  cm)  did  not 
dominate  the  balance,  rather  It  was  the  Intermediate  thicknesses  of 
young  Ice  (10-80  cm)  that  were  responsible  for  over  half  the  total  loss. 

Total  Ice  production  within  the  region  <f>  was  also  dramatically 
larger  when  contributions  from  young  Ice  were  considered.  Figure  29a 
shows  annual  Ice  production  In  each  of  the  thickness  categories.  The 
mass  produced  In  the  0-10  cm  category,  for  example,  was  equivalent  to  a 
15.4  cm  thick  layer  covering  the  entire  region.  For  the  year  <f> 
totaled  about  89  cm,  In  contrast  to  the  zero  mass  balance  expected  for 
3  m  Ice.  Growth  In  areas  of  first-year  Ice  accounted  for  essentially 
all  of  the  net  Ice  production.  Despite  the  extremely  large  growth  rates 
In  new  leads,  only  about  15%  of  the  first-year  Ice  production  occurred 
In  the  0-10  category,  reflecting  the  very  small  area  occupied  by  open 
water  during  the  winter.  As  expected,  the  greatest  production  took 
place  following  the  fall  freezeup.  Seasonal  variations  In  the  flux  of 
salt  to  the  ocean  paralleled  the  Ice  production.  Net  salt  Input  began 
In  mid- September ,  the  observed  time  of  minimum  salt  content  In  the  mixed 
layer.  Without  young  Ice,  the  salt  content  of  the  mixed  layer  would  not 
begin  to  Increase  significantly  until  early  December. 

Numerical  experiments  were  also  conducted  to  determine  how 
uncertainties  In  specified  strains  might  affect  regional  heat  and  mass 
fluxes .  It  was  found  that  the  average  divergence  of  the  Ice  was  poorly 
correlated  with  monthly  or  annual  heat  flux  totals,  but  that  <f>  was 
related  to  the  time  average  of  the  strain  components.  Ice  typically 
alternates  between  diverging  and  converging  motions  over  periods  of 
days,  producing  continual  changes  In  the  amount  of  open  water.  Ice 
production,  and  ridging  activity.  It  appears  that  these  short-term 
variations  In  strain  are  more  Important  to  the  thermodynamics  than  the 


longer  term  averages.  The  tendency  for  <£>  to  Increase  with  Increased 
divergence  suggests  that  variations  In  the  export  of  Ice  from  the  Arctic 
Basin  may  not  have  a  great  Impact  on  total  Ice  volume  because  of  the 
compensating  changes  In  Ice  production. 

These  calculations  have  shown  that  thin  Ice  and  open  water 
resulting  from  differential  Ice  movement  cause  Interactions  between  the 
ocean  and  atmosphere  to  be  much  more  vigorous  than  would  be  the  case  for 
a  uniform,  static  Ice  cover.  Further,  It  Is  evident  that  care  must  be 
exercised  when  Interpreting  local  measurements  made  In  areas  of 
nonuniform  Ice.  Failure  to  consider  thickness  variations  and  how  they 
affect  area  averages  can  lead  to  substantial  errors  In  understanding 
large-scale  Interactions  between  the  Ice,  ocean,  and  atmosphere. 
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